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WELL-POSEDNESS FOR THE MOTION OF AN 
INCOMPRESSIBLE LIQUID WITH FREE SURFACE BOUNDARY 

o 
o 
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[j^ , Abstract. We study the motion of an incompressible perfect liquid body in vacuum. This can be thought 

of as a model for the motion of the ocean or a star. The free surface moves with the velocity of the liquid 
and the pressure vanishes on the free surface. This leads to a free boundary problem for Euler's equations, 
where the regularity of the boundary enters to highest order. We prove local existence in Sobolev spaces 
assuming a "physical condition", related to the fact that the pressure of a fluid has to be positive. 

< 

! 1. Introduction 

We consider Euler's equations describing the motion of a perfect incompressible fluid in vacuum: 

(1.1) {dt + V^dk)vj+djP = 0, j = l,...,n in V, 

m _ 

^ ■ where di = d/dx^ and V = Uo<t<T {t} x Vt, Vt C M". Here V''= = Vk, and we use the convention 

Tjlj- ! that repeated upper and lower indices are summed over. Vis the velocity vector field of the fluid, p is 

CD ' the pressure and T>t is the domain the fluid occupies at time t. We also require boundary conditions on 

r~| ■ the free boundary dV = Uo<t<T {t} x dVt; 

(1.3) p = 0, on dV, 

(1.4) {dt + V''dk)\aveT{dV), 

\^ • Condition (1.3) says that the pressure p vanishes outside the domain and condition (1.4) says that the 

^ boundary moves with the velocity V of the fluid particles at the boundary. 

Given a domain Dq C M", that is homeomorphic to the unit ball, and initial data vq, satisfying 
the constraint (1.2), we want to find a set 2? = Uo<t<r {t} x 'Dt, Vt C and a vector field v solving 
(1.1)-(1.4) with initial conditions 

(1.5) {x; (0,x) e P} = Po, and v = vq, on {0} x 

Let M be the exterior unit normal to the free surface dDf Christodoulou[C2] conjectured that the 
initial value problem (1.1)-(1.5), is well posed in Sobolev spaces if 

(1.6) VArp<-co<0, on dV, where Vj^f = N'd^i. 
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Condition (1.6) is a natural physical condition since the pressure p has to be positive in the interior 
of the fluid. It is essential for the well posedness in Sobolev spaces. A condition related to Rayleigh- 
Taylor instability in [BHL,W1] turns out to be equivalent to (1.6), see [W2]. Taking the divergence of 
(1.1) gives: 

(1.7) -Ap={djV'')dkV^, in Vt, p = 0, on dVt 

In the irrotational case, when curlvij = diVj — djVi = 0, then Ap < so p > and (1.6) holds by the 
strong maximum principle. Furthermore Ebin [El] showed that the equations are ill posed when (1.6) 
is not satisfied and the pressure is negative and Ebin [E2] announced an existence result when one adds 
surface tension to the boundary condition which has a regularizing effect so (1.6) is not needed then. 

The incompressible perfect fluid is to be thought of as an idealization of a liquid. For small bodies 
like water drops surface tension should help holding it together and for larger denser bodies like stars 
its own gravity should play a role. Here we neglect the influence of such forces. Instead it is the 
incompressibility condition that prevents the body from expanding and it is the fact that the pressure 
is positive that prevents the body from breaking up in the interior. Let us also point out that, from 
a physical point of view one can alternatively think of the pressure as being a small positive constant 
on the boundary instead of vanishing. What makes this problem difficult is that the regularity of the 
boundary enters to highest order. Roughly speaking, the velocity tells the boundary where to move 
and the boundary is the zero set of the pressure that determines the acceleration. 

In general it is possible to prove local existence for analytic data for the free interface between two 
fluids. However, this type of problem might be subject to instability in Sobolev norms, in particular 
Rayleigh-Taylor instability, which occurs when a heavier fluid is on top of a lighter fluid. Condition 
(1.6) prevents Rayleigh-Taylor instability from occurring. Indeed, if condition (1.6) is violated Rayleigh- 
Taylor instability occurs in a linearized analysis. 

Some existence results in Sobolev spaces were known in the irrotational case, for the closely related 
water wave problem which describes the motion of the surface of the ocean under the influence of earth's 
gravity. The gravitational field can be considered as uniform and it reduces to our problem by going 
to an accelerated frame. The domain Vt is unbounded for the water wave problem coinciding with a 
half-space in the case of still water. Nalimov[Na] and Yosihara[Y] proved local existence in Sobolev 
spaces in two space dimensions for initial conditions sufficiently close to still water. Beale, Hou and 
Lowengrab[BHL] have given an argument to show that problem is linearly well posed in a weak sense 
in Sobolev spaces, assuming a condition, which can be shown to be equivalent to (1.6). The condition 
(1.6) prevents the Rayleigh-Taylor instability from occurring when the water wave turns over. Finally 
Wu[Wl,2] proved local existence in the general irrotational case in two and three dimensions for the 
water wave problem. The methods of proofs in these papers uses that the vector field is irrotational to 
reduce to equations on the boundary and do not generalize to deal with the case of nonvanishing curl. 

We consider the general case of nonvanishing curl. With Christodoulou [CL] we proved local a 
priori bounds in Sobolev spaces in the general case of non vanishing curl, assuming (1.6) hold initially. 

Usually if one has a priori estimates, existence follows from similar estimates for some regularization 
or iteration scheme for the equation, but the sharp estimates in [CL] use all the symmetries of the 
equations and so only hold for perturbations of the equations that preserve the symmetries. In [LI] 
we proved existence for the linearized equations, but the estimates for the solution of the linearized 
equations looses regularity compared to the solution we linearize around, so existence for the nonlinear 
problem does not follow directly. Here we use improvements of the estimates in [LI] together with the 
Nash-Moser technique to show local existence for the nonlinear problem in the smooth class: 



Theorem 1.1. Suppose that vq and BVq in (1-5) are smooth, Vq is diffeomorphic to the unit ball, and 
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that (1.6) hold initially when t = 0. Then there is a T > such that (1.1)-(1.5) has a smooth solution 
for <t <T, and (1-6) hold with cq replaced by co/2 for <t <T. 

In [CL] we proved local energy bounds in Sobolev spaces. It now follows from the bounds there that 
the solution remains smooth as long as it is and the physical condition (1.6) hold. The existence 
for smooth data now implies existence in the Sobolev spaces we considered in [CL]. Moreover, method 
here also works for the compressible case [L2,L3]. 

Let us now describe the main ideas and difficulties in the proof. In order to construct an iteration 
scheme wc must first introduce some parametrization in which the moving domain becomes fixed and 
express Eulcr's equations in this fixed domain. This is achieved by the Lagrangian coordinates given 
by following the flow lines of the velocity vector field of the fluid particles. 

In [LI] we studied the linearized equations of Euler's equations expressed in Lagrangian coordinates. 
We proved that the linearized operator is invertible at a solution of Euler's equations. The linearized 
equations become an evolution equation for what wc called the normal operator, (2.17). The normal 
operator is unbounded and not elliptic but it is symmetric and positive on divergence free vector fields if 
(1.6) hold. This leads to energy bounds and existence for the linearized equations follows from a delicate 
regularization argument. The solution of the linearized equations however looses regularity compared 
to the solution we linearize around so existence for the nonlinear problem docs not follow directly from 
an inverse function theorem in a Banach space but wc must use the Nash-Moscr technique. 

We first define a nonlinear functional whose zero will be a solution of Euler's equations expressed 
in the Lagrangian coordinates. Instead of defining our map by the left hand side of (1.1) and (1.2) 
expressed in the Lagrangian coordinates we let our map be given by the left hand side of (1.1) and 
we let pressure be implicitly defined by (1.7) satisfying the boundary condition (1.3). This is because 
one has to make sure that the pressure vanishes on the boundary at each step of an iteration or else 
the linearized operator is ill posed. One can see this by looking at the irrotational case where one 
gets an evolution equation on the boundary. If the pressure vanishes on the boundary then one has an 
evolution equation for a positive elliptic operator but if it does not vanish on the boundary there will 
also be some tangential derivative, no matter how small coefficients they come with the equation will 
have exponentially growing Fourier modes. 

In order to use the Nash-Moser technique one has to be able to invert the linearized operator in a 
neighborhood of a solution of Euler's equations or at least do so up to a quadratic error [Ha]. In this 
paper we generalize the existence in [LI] so the linearized operator is invertible in a neighborhood of a 
solution of Euler's equations and outside the class of divergence free vector fields. This does present a 
difficulty because the normal operator, introduced in [LI] is only symmetric on divergence free vector 
fields and in general it looses regularity. Overcoming this difficulty requires two new observations. The 
first is that, also for the linearized equations there is an identity for the curl that gives a bound that is 
better than expected. The second is that one can bound any first order derivative of a vector field by 
the curl, the divergence and the normal operator times one over the constant cq in (1.6). Although the 
normal operator is not elliptic on general vector fields it is elliptic on irrotational divergence free vector 
fields and in general one can invert it if one also have bounds for the curl and the divergence. 

The methods here and in [CL] arc on a technical level very different but there are philosophical 
similarities. First wc fix the boundary by introducing Lagrangian coordinates. Secondly, wc take the 
geometry of the boundary into account. Here in terms of the normal operator and Lie derivatives with 
respect to tangential vector fields and in [CL] in terms of the second fundamental form of the boundary 
and tangential components of the tensor of higher order derivatives. Thirdly, we use interior estimates 
to pick up the curl and the divergence. Lastly, we get rid of a difficult term, the highest order derivative 
of the pressure, by projecting. Here we use the orthogonal projection onto divergence free vector fields 
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whereas in [CL] we used the local projection of a tensor onto the tangent space of the boundary. 

The paper is organized as follows. In section 2 we reformulate the problem in the Lagrangian 
coordinates and give the nonlinear functional which a solution of Euler's equations is a zero of and we 
derive the linearized equations in this formulation. In section 2 we also give an outline of the proof 
and state the main steps that we will prove. The main part of the paper, sections 3 to 13 are devoted 

to proving existence and tame energy estimates for the inverse of the linearized operator. Once this is 
proven, the remaining sections 14 to 18 are devoted to setting up the Nash-Moser theorem we are using. 



2. Lagrangian coordinates and the linearized operator 



Let us first introduce the Lagrangian coordinates in which the boundary becomes fixed. By a 
scaling we may assume that Vq has the volume of the unit ball CI and since we assumed that Vq is 
diffeomorphic to the unit ball we can, by a theorem in [DM], find a volume preserving diffeomorphism 
fo-.Q—y Vq, i.e. det (dfo/dy) = 1. Assume that v{t,x), p{t,x), {t,x) £ V are given satisfying the 
boundary conditions (1.3)-(1.4). The Lagrangian coordinates x = x{t,y) = ft{y) are given by solving 

(2.1) '^^^ = V{t,x{t,y)), x{0,y) = My), yen 

Then : -^^ 2?^ is a volume preserving diffeomorphism, if divF = 0, and the boundary becomes fixed 
in the new y coordinates. Let us introduce the material derivative: 



(2.2) A = I 



_ d_ 

i=constant dt 



+ V , 

x=constant 



The partial derivatives di = d/dx^ can then be expressed in terms of partial derivatives da = d/dy"- 
in the Lagrangian coordinates. We will use letters a,b,c,...,f to denote partial differentiation in the 
Lagrangian coordinates and i,j,k, ... to denote partial differentiation in the Eulerian frame. 

In these coordinates Euler's equation (1.1) become 

(2.3) D^x, + dip = 0, (t, y) G [0, T] x n, 

where now Xi = Xi(t, y) and p = p(t, y) are functions on [0, T] x Vt, Dt is just the partial derivative with 
respect to t and di = {dy"" / dx^)da, where da is differentiation with respect to y". (1.7) become 

(2.4) Ap+{diV'')dkV' = 0, p =0, where V' = Dtx\ 



" dx^ dx^ 



Here 

(2.5) Ap = Y^^dfp=K~^da{Kg''^dbp) where 9ab = ^zj Q^a Qyb ■ 

g°'^ is the inverse of the metric gab and k = det (dx/dy) = V det g. The initial conditions (1.5) become 

(2.6) xl^^Q = /o, = vo 



Christodoulou's physical condition (1.6) become 

(2.7) V^rp<-co<0, on dn, where VAf=M'd^,. 
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This is needed in the proof for the normal operator (2.17) to be positive which leads to energy bounds. 
In addition to (2.7) we also need to assume a coordinate condition having to do with that we are 
looking for a solution in the Lagrangian coordinates and we are starting by composing with a particular 
diffeomorphism. The coordinate condition is 

71 

(2.8) \dx/dyf + \dy/dxf < c?, ^ (|5«^| + \9ab\) < ncj, 

a, 6=1 

where \dx/dy\^ = Yll a=i^^^^ / ^V^T ■ This is needed for (2.5) to be invertible. We note that the second 
condition in (2.8) follows from the first and the first follows from the second with a larger constant. 
We remark that this condition is fulfilled initially since we are composing with a diffeomorphism. 
Furthermore, for solution of Euler's equations, divF = 0, so the volume form k is preserved and hence 
an upper bound for the metric also implies a lower bounded for the eigenvalues and an upper bound 
for the inverse of the metric follows. However, in the iteration, we will go outside the divergence free 
class and hence we must make sure that both (2.7) and (2.8) hold at each step of the iteration. We will 
prove the following theorem: 

Theorem 2.1. Suppose that initial data (2.6) are smooth, vq satisfy the constraint (1-2), and that (2.7) 
and (2.8) hold whent = 0. Then there isT > such that (2.3)-(2.4) has a solution x,p £ C°° {[0,T]xTl). 
Furthermore, (2.7)-(2.8) hold, for <t <T, with cq replaced by co/2 and c\ replaced by 2c\. 

Theorem 1.1 follows from Theorem 2.1. In fact, the assumption that is diffeomorphic to the unit 
ball, together with that one then can find a volume preserving diffeomorphism guarantees that (2.8) 
hold initially. Once, we obtained a solution to (2.3)-(2.4), we can hence follow the flow lines of V in 
(2.1) and this defines a diffeomorphism of [0,T] x Q to P, so we obtain smoothness of F as a function 
of (t, x) from the smoothness function of (t, y). 

In this section we first define a nonlinear functional whose zero is a solution of Euler's equations, 

(2.9) -(2.13). Then we derive the linearized operator in Lemma 2.2. The existence will follow from the 
Nash-Moser inverse function theorem, once we proven that the linearized operator is invertible and so 
called tame estimates for the inverse stated in Theorem 2.3. Proving that the linearized operator is 
invertible away from a solution of Euler's equations and outside the divergence free class is the main 
difficulty of the paper. This is because the normal operator (2.17) is only symmetric and positive 
within the divergence free class and in general it looses regularity. In order to prove that the linearized 
operator is invertible and estimates for its inverse we introduce a modification (2.31) of the linearized 
operator that preserves the divergence free condition, and first prove that the modification is invertible 
and estimates for its inverse, stated in Theorem 2.4. The difference between the linearized operator 
and the modification is lower order and the estimates for the inverse of the modified linearized operator 
leads to existence and estimates also for the inverse of the linearized operator. 

Proving the estimates for the inverse of the modified linearized operator, stated in Theorem 2.4, 
takes up most of the paper, sections 3 to 13. In this section we also derive certain identities for the 
curl and the divergence, see (2.29)-(2.30), need for the proof of Theorem 2.4. Here we also transform 
the vector field to the Lagrangian frame and express the operators and identities in the Lagrangian 
frame, see Lemma 2.5. The estimates in Theorem 2.4 will be derived in the Lagrangian frame since 
commutators of the normal operator with certain differential operators are better behaved in this frame. 

In section 3, we introduce the orthogonal projection onto divergence free vector field and decompose 
the modified linearized equation into a divergence free part and an equation for the divergence. This is 
needed to prove Theorem 2.4 because the normal operator is only symmetric on divergence free vector 
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fields and in general it looses regularity. However, we have a better equation for the divergence which 
will allow us to obtain the same space regularity for the divergence as for the vector field itself. 

In section 4 we introduce the tangential vector fields and Lie derivatives and calculate commutators 
between these and the operators that occur in the modified linearized equation, in particular the normal 
operator. In section 5 we show that any derivative of a vector field can be estimated by derivatives of the 
curl and of the divergence, and tangential derivatives or tangential derivatives of the normal operator. 
In section 6 introduce the L°° norms that we will use and state the interpolation inequalities that we will 
use. In section 7 and 8 we give the tame L^oo and L°° estimates for the Dirichlet problem. In section 
9 we give the equations and estimates for the curl that we will use. In section 10 we show existence for 
the modified linearized equations in the divergence class. In section 11 we give the improved estimates 
for the inverse of the modified linearized operator within the divergence free class. These are needed in 
section 12 to prove existence and estimates for the inverse of the modified linearized operator. Finally 
in section 13 we use this to prove existence and estimates for the inverse of the linearized operator. 

In section 14 we explain what is needed to ensure that the physical and coordinate conditions (2.7) 
and (2.8) continue to hold. In section 15 we summarize the tame estimates for the inverse of the 
linearized operator in the formulation that will use with the Nash-Moser theorem. In section 16 we 
derive the tame estimates for the second variational derivative. In section 17 we give the smoothing 
operators needed for the proof of the Nash-Moser theorem on a bounded domain. Finally, in section 18 
we state and prove the Nash-Moser theorem in the form that we will use. 

Let us now define the nonlinear map, that we will use to find a solution of Euler's equations. Let 

(2.9) ^i = ^i{x) = Dlxi + dip, where di = idyydx')da, 
and p = *(x) is given by solving 

(2.10) Ap = -{diV'')dkV\ p\9a=0, where V = Dtx. 
A solution to Euler's equations is given by 

(2.11) $(x) = 0, for 0<t<T, xl^^ = fo, DtX^^^ = vq 

We will find T > and a smooth function x satisfying (2.11) using the Nash-Moser iteration scheme. 

First we turn (2.11) into a problem with vanishing initial data and a small inhomogeneous term 
using a trick from [Ha] as follows. It is easy, to construct a formal power series solution xq as t ^ 0: 

(2.12) D^^xo)l^^ = 0, k>0, xol^^ = fo, Axo|,^o = «o 

In fact, the equation (2.10) for the pressure p only depends on one time derivative of the coordinate x 
so commuting through time derivatives in (2.10) gives a Dirichlet problem for D^p depending only on 
D^x, for m < k + 1 and Dfp, for i < k — 1. Similarly commuting through time derivatives in Euler's 
equation, (2.11), gives D^'^^x in terms of for m < k, and D^p, for £ < k. We can hence construct 

a formal power series solution in i at i = and by a standard trick we can find a smooth function xq 
having this as it power series, see section 10. We will now solve for u in 

(2.13) ^u) = ^iu + xo)-^xo) = Fs-Fo = fs, u^^^ = Dtu^^^ = 

where Fs is constructed as follows. Let Fq = $(.xo) and let Fs{t,y) = Fo{t — 5,y), when t > S and 
Fs{t,y) = 0, when t < 6. Then Fs is smooth and fs = Fs — Fo tends to in C°° when 5^0. 
Furthermore, fs vanish to infinite order as t ^ 0. Now, ^(0) = so it will follow from the Nash- 
Moser inverse function theorem that i(u) = fs has a smooth solution -u if 5 is sufficiently small. Then 
X = u + xo satisfies (2.11) for <t < S. 

In order to solve (2.11) or (2.13) we must show that the linearized operator is invertible. Let us 
therefore first calculate the linearized equations. Let 5 be the Lagrangian variation, i.e. derivative w.r.t. 
some parameter r when {t, y) are fixed. We have: 
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Lemma 2.2. Let x = x{r, t, y) he a smooth function of (r, t, y) ^ K = [— e, e\ x [0, T] x Vl, e > such 
that x\^_Q = X. Then ^{x) is a smooth function of (r,t,y) G K, such that d^{x)/dr\^_^ = ^'{x)dx, 
where 6x = dx/dr\^_^ and the linear map Lq = <^'{x) is given by 

(2.14) ^'{x)5xi = D^Sxi + {dkdip)5x^ + diSpo + di [Spi - Sx^dkp) , 
where p satisfies (2.10) and Spi, i = 0,1, are given by solving 

(2.15) A{6pi-5x''dkp) =0, ^Pi\dn=^^ 

(2.16) A5po = -2{dkV')di{5V'' - Sx'diV'^), Spo\g^= 0, 

where Sv = DfSx. Here, the normal operator 

(2.17) ASxi = -di{dkP Sx'' - dpi) 

restricted to divergence free vector fields is symmetric and positive, in the inner product {u, w) = 
Jj,^ dx, if the physical condition (2.7) hold. 

Proof. That $(x) is a smooth function follows from that the solution of (2.10) is a smooth function if 
X is, see section 16. Let us now calculate $'(x). Since [S,d/dy°'] = it follows that 

(2.18) [S^9,] = {^%)^a-i9^Sx^)9u 

where we used the formula for the derivative of the inverse of a matrix = -A-^{5A)A-^. It 

follows that [d—Sx'^di,di] = {5—5x^di is the Eulerian variation ). Hence 

(2.19) (^$, - Sx^dk^^ = DjSxi - {dkD^Xi)6x^ + ^^ {5p - Sx'^dtp) , where 

(2.20) A{6p - dx'^dkp) = {5- 5x^dk)Ap = -2{dkV')di{6V'' - Sx'diV''), 6p\g^= 0. 

The symmetry and positivity of A were proven in [LI], see also section 3 here. □ 

In order to use the Nash-Moser iteration scheme to obtain a solution of (2.13) we must show that 
linearized operator is invertible and that the inverse satisfies tame estimates: 

Theorem 2.3. Let 

(2.21) lll«ll|a,fc= sup \\u{t,-)\\a,oo + ■■■ + \\D'lu{t,-)\\a,o^ 

0<t<T 

where ||u(t, ■)||a,oo o,f^ the Holder norms in see (17.1). 

Suppose that (2.7) and (2.8) hold initially, where p is given by (2.10), and let xq G C°°([0,r] X O) 
satisfy (2.12). Then there is a Tq = T(xo) > 0, depending only on upper bounds for |||a:;o|||4,2; Cg and 
c\, such that the following hold. If x E C°°{jO,T] x p is defined by (2.10), 

(2.22) T<ro, |||x-xo||k2 < 1, and {x - xo)\^^^ = Dt{x - xo)\^^^ = 

then (2.7) and (2.8) hold for < t <T with cq replaced by co/2 and ci replaced by 2c\. Furthermore, 
linearized equations 

(2.23) ^'{x)Sx = 5$, in [0, T] x Q, 6x\^^^ = Dt 5x\^^^ = 0. 
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where G C°°([0,T] x $7) has a solution Sx G C°°([0,T] x . T/ie solution satisfies the estimates 



(2.24) \\\Sx\\\a,2 < Ca{\m\\\a+ro+2,0 + || |<5<^ || 1 1,0 |||x " Xq || U+ro+6,2) , « > 

where Ca = Ca(xo) is bounded when a is bounded, and in fact depends only on upper bounds for 
|||xo|||o+ro+6,2; and ci. Here ro = [w/2] + 1, where n is the number of space dimensions. 

Furthermore $ is twice differentiable and the second derivative satisfies the estimates 
(2.25) 

|||$"(x)((5x,ex)|||a,o < C'a(|||<5a;|||a+4,i|||ex|||2,i + |||<5a;|||2,i|||ex|||a+4,i + |||a;-xo|||a+4,i|||5a;|||2,i|||ex|||2,i) 



The proof of Theorem 2.1 follows from Theorem 2.3 and Proposition 18.1. In Theorem 2.3 
we use norms that only has two time derivatives and our Nash-Moser theorem, Proposition 18.1, 
gives a solution of (2.13) u G C^([0, T], C°°(i7)). However, additional regularity in time follows 
from differentiating the equations with respect to time. In fact, if a; G C'^ ([0, T], C°°(0)) then 
D^x = —dip G C'^~^([0,T],C°°(J7)), since (2.10) only depends on one time derivative of x, see the 
proof of Lemma 6.7, and it follows that x G ^=+^([0, T], C°°{n)) . 

Theorem 2.3 follows from Lemma 14.1, Proposition 15.1 and Proposition 16.1. The main point 
being existence for (2.23) and the tame estimate (2.24) given in Proposition 15.1. We will now discuss 
how to prove existence and estimates for the linearized equations. The terms {dkdip)5x^ and diSpo 
in (2.14) are order zero in 5x and DfSx. The last term is a positive symmetric operator but only on 
divergence free vector fields and in general it is an unbounded operator that looses regularity. In general 
6x is not going to be divergence free but we will derive evolution equations for the divergence and the 
curl of Sx, that gain regularity. These evolution equations comes from that the divergence and the 
curl of the velocity v are conserved expressed in the Lagrangian coordinates for a solution of Euler's 
equations, $(x) = 0. In fact, since [Dt,di] = —{diV^)dk it follows from (2.9) that 

(2.26) DtdivV = div$, Cd^ curl?; = curl$ 

where cuvlvij =diVj— djVi and is the space time Lie derivative with respect to Dt = (l,V): 

(2.27) £D,<Tij = Dt aij + idiV^)aij + {djV^)aa 

restricted to the space components. Expressing the two form a in the Lagrangian frame this is just the 
time derivative: 

(2.28) Dt{aiaiaij) =aiaiCD,c7ij, where < = SxVSy" 

We have the following evolution equations for the divergence and the curl of the linearized operator 

(2.29) dw{^\x)Sx) = D^divSx + {diSx^)dk^\ 

(2.30) curl($'(x)(5x) = Co, curl (A Sx - Sx'^dvk) + {diSx'')dj^k - {djSx^)di^k 

In fact, since [5, 9j] = —{diSx'^)dk and [Dt,di\ = —{diV^)dk it follows that SdivDtx = DtdivSx so 
by (2.26) D'^dYvSx = 5div$ and (2.29) follows. To prove (2.30) we note that [S,aiaidi] = [S,aida] = 
{Sal)da = {dbSx^)da = ala^{dkSx^)di so 

(2.31) S[alaicwclvij) = alal^cuilSvij + {djSx^)diVk - {diSx^)dkVk) = aj^a^ curl ((5 v - SxkdV'')ij 
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where cutI{Sv— 5xkdV'')ij = di{5vj—5x^djVk)—dj{5vi—5x^diVk) and (2.30) follows since by (2.26)-(2.28) 



(2.32) Cd, curl {5v - Sx^dV'') = curl (5$ - Sx^d^''). 

In [LI] we proved existence and estimates for the inverse of tlie linearized operator at a solution 
of Euler's equations and within the divergence free class. We only inverted ^'{x)5x = 5$ when 5$ 
was divergence free and <t{x) = 0, in which case by (2.29) Sx is also divergence free. In order to use 
the Nash-Moscr iteration scheme we will show that the linearized operator is invcrtiblc away from a 
solution of Euler's equations and outside the divergence free class. This does present a problem since 
the normal operator is only symmetric on divergence free vector fields so for general vector fields we 
loose a derivative. In order to recover this loss we will use that one has better evolution equations for 
the divergence and for the curl that do not loose regularity. (2.29)-(2.30), says that wc can get bounds 
for the divergence and the curl of dx if we have bounds for all first order derivatives of 5x. In fact 
(2.29)-(2.30) can be integrated even without knowing a bound for first order derivatives of Df Sx. 

We will now first modify the linearized operator so as to remove the term {didx'')dk^'' in (2.29) 
without making (2.30) worse. (2.29) without this term will give us an evolution equation that allows us 
to control the divergence. This together with that the normal operator (2.17) is symmetric and positive 
on divergence free vector fields will give us existence for the inverse of the modified linearized operator. 
The modified linearized operator is given by 

(2.33) LiSx' = ^'{x)6x' - Sx'^dk^' + Sx'div^ 

= D^6x^ - {dkD^x')6x'' + di {6pi - Sx'^dkp) + 5x' div$ + diSpo 

It follows from (2.29) that 

(2.34) div (Life) = D"^ divfe + div$ divfe 

The operator Li reduces to the linearized operator Lq = $'(x) when $(x) = and the difference Li — Lq 
is lower order. Furthermore, Li preserves the divergence free condition. We will first prove existence for 
the inverse of the modified linearized operator and the existence of the inverse of the linearized operator 
follows since the difference it is a lower order. The main part of the manuscript is devoted to proving 
the following existence and energy estimates: 

Theorem 2.4. Suppose that x is smooth and that the physical condition (2.7) and the coordinate 
condition (2.8) hold forO<t<T. Then 

(2.35) LiSx = (5$, < t < r, Sx\^^^ = Dt Sx\^^^^ = 

has a smooth solution 6x if 6^ is smooth. 

Furthermore, there are constants K4 depending only on upper bounds for T, Cq^ , ci, r and |||a;|||4_2 
such that the following estimates hold, //div (5$ = then divSx = and 

(2.36) \\Dtdx\\r + \\Sx\\r < K4 [ {\\6^r + \\\x\\\r+3,i\\S^o) dr, r>0. 

Jo 

If div 6^ = 0, curl(5# = and = then 

(2.37) + ||Dtfe||r + \\Sx\\r + Co\\6x\\r+l 

<Ki! (||A5$||r + ||'5$||r + |||a;|||r+3,2(||A5$||o + ||'5$||o))cir, r>0 
^0 
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In general 



(2.38) II Afe||^_i + ||(5x||^ < ^4 / (||5$||^ + |||ar|||^+3,2||<5$||i)rfT, r>l 

Jo 

Here |||aj|||r,fc is as in Theorem 2.3 and 

(2.39) ||fe||. = \\6x{t,-)\\r = E ( / \dySx{t,y)\'dyy 



1/2 



|a| <r 



The proof of the existence for (2.23) and the tame estimate (2.24) for the inverse of the hnearized 
operator in Theorem 2.3 follows from Theorem 2.4. In fact, since the difference (Li — ^'{x))Sx = 0{5x) 
is lower order, the estimate (2.38) will then allow us to get existence and the same estimate also for 
the inverse of the linearized operator (2.23), by iteration. In (2.38) we only have estimates for one time 
derivative, but we get estimates for an additional time derivative from also using the equation. The 
estimates for (2.23) so obtained then gives the L°° estimates (2.24) by also using Sobolev's lemma. 

The proof of Theorem 2.4 takes up most of the manuscript. The proof (2.36) uses the symmetry and 
positivity of the normal operator (2.17) within the divergence free class. This leads to energy estimates 
within the divergence free class. The proof of (2.37) is obtained by first differentiating the equation 
with respect to time and then using that a bound for two time derivatives also gives a bound for the 
normal operator (2.17) using the equation. The normal operator is not elliptic acting on general vector 
fields. However, it is elliptic acting on divergence and curl free vector fields and in general one can 
invert it and gain a space derivative if one also has bounds for the curl and the divergence, see Lemma 
5.4. Here we also need to use the improved estimate for the curl coming from (2.30). To prove (2.38) we 
first subtract of a vector field picking up the divergence. The equation for the divergence from (2.34): 

(2.40) divSx + div$ divSx = div5$ 

is just an ordinary differential equation that do not loose regularity and in fact the estimates for (2.40) 
gain an extra time derivative compared to the estimate (2.36). Once we control the divergence we use 
the orthogonal projection onto divergence free vector fields to obtain an equation for the divergence free 
part by projecting the equation (2.35), see section 3. The equation so obtained is of the form (2.35) 
with divS^ = and 6^ depending also on the divergence div6x that we just calculated. The interaction 
term coming from the divergence part looses a space derivative but it is in the form of a gradient so we 
can recover this loss by using the gain of a space derivative in (2.37). 

In order to prove the energy estimates needed to prove Theorem 2.4 one has to express the vector 
fields in the Lagrangian frame, see (2.43). Theorem 2.4, expressed in the Lagrangian frame, follows 
from Theorem 10.1, Theorem 11.1 and Theorem 12.1. Below, we will express the equation (2.35) in the 
Lagrangian frame and in section 3 we outline the main ideas of how to decompose the equation into a 

divergence free part and an equation for the divergence using the orthogonal projection onto divergence 
free vector fields and we show the basic energy estimate within the divergence free class. 

As described above we now want to invert the modified linearized operator (2.35) by decomposing 
it into an operator on the divergence free part and the ordinary differential equation (2.40) for the 

divergence. Hence we first want to be able to invert Li in the divergence free class. The normal operator 
A, the third term on the second row in (2.33), maps divergence free vector fields onto divergence free 
vector fields. We also want to modify the time derivative by adding a lower order term so it preserves 

10 



the divergence free condition. Let the Lie derivative and modified Lie derivative with respected to the 
time derivative acting on vector fields be defined by 

(2.41) CdM = DtSx' - {dkV')5x^, and CdM = J^dM + divF fe* 
As before, is the space time Lie derivative restricted to the space components. Then 

(2.42) div jO-DtSx = Dt divSx, where Dt = Dt + divV 

i.e. Dtf = Dtf + {dwV)f. 

This is easier to see if we express the vector field in the Lagrangian frame. Let 

(2.43) = -^Sx^ 
Then, 

(2.44) Dt Sx' = Dt {W''dxydy'') = {DtW^)dx^ / dy^ + W^dV^/dy^ = {DtW^)dx' / dy'' + 5x''dkV' 
and multiplying with the inverse dy°' jdx^ gives 

(2.45) DtW'' = ^-^^LDM\ and DtW =^-£-LdM ■ 
With K = det (dx/dy), we have 

(2.46) W^" = D^W'' = DtW + (divF)W^" = K-^DtinW) 
since DtK = KdivV, see [LI]. Since the divergence is invariant 

(2.47) divfe = divVF = K-^daiKW") 
it therefore follows that 

(2.48) divAw^ = AdivW^ 

The idea is now to replace the time derivatives Dt in (2.33) by C^t or equivalcntly express Li in 
the Lagrangian frame and use the modified time derivatives Dt. Expressing the operator Li in the 
Lagrangian frame we get: 

Lemma 2.5. Let W = DtW and W = D^^W . Then we can write (2.35) as L^W = F, where W is 
given by (2.43), F« = ^^dy°-/dx^ and 

(2.49) L^W" = W'' + - B{W,WY, B{W,Wf = BQW'' + B{W''. 
Here 

(2.50) QabAW'' = -da {{d.p)^ - gi) , dxvAW = 

(2.51) QabB^W'' = a{Dtgac-uJac-CTgac)W^-daq^, div B^W = -a^diwW 

(2.52) gahBrW'' = -{Dtgac - ^ac - 2&gac)W'' - 5„g2, div^iW^ = 2adivW, 

11 



where qi, for i = 1, 2, 3 are given by solving the Dirichlet problem qilg^ = where Aqi are given by the 
equations for the divergences above, cr = In a = D^a = divF, a = D^a and 

(2.53) A 9ah=Q:^^ {diV^ + d^Vi) , a;„6 = — ^ {d^Vj - djVi) . 
We have 

(2.54) div(LiW^) = D'^diYW + adiYW 

Let L]Wa = QabLiW^, Wa = QahW^ and Wa = Wa- {iVab + c^9ab)W'. Then 

(2.55) curl (InW) = Dt curly} + curlS4W" 

(2.56) curl (LiW) = Dt curM + curlSg + curlSgW" 

where B^Wa = {DtLUab + (^9ab)W^ , B^Wa = -{uJab + c^gabW^ andB^Wa = -CT{Dt gab-i^ab- (^9ab)W' . 
Furthermore Lq = #' (x) expressed in the Lagrangian frame is given by 

(2.57) LoW = LiW" - BsW, where BsW = -W^Vc^'' + W div^ 

where Vc is covariant differentiation with respect to the metric gab and = ^^dy"' /dx^ , i.e. Vc^"^ = 
{dx^ldy''){dy<'/dxi)di^^. 

Proof. Differentiating (2.44) once more gives 

(2.58) Dl5x' - {dkDtV')8x^ = {DlW^)dx' /dy^ + 2{DtW'')dV' / dy^ 
It follows that 

(2.59) 5y« ^ * V t y ; dy"' dy^ * ^ y t ' Qyb Qya » J 

= gabDlW' + (A 9ab - C^ab) AW^*" 

It follows from (2.33) that 

(2.60) gabLiW'' = gabD^W'' - {{deP)W^ -q) + (A5ac - u;ac)DtW'^ + agabW' 

= Dt{gabDtW'^ - LOabW'') - da{{d,p)W- - q) + DtUabW' + agabW' 

where q = Sp is chosen so that the divergence is equal to divLiVF = divVF + div A divF in order 
for it to be consistent with (2.34). We have = {Dt + divy)(A + divF) = D^ + 2a Dt + a'^ + '& = 
D^ + 2&Dt + a-&^ so 

(2.61) D^ = - 2abt + cr'^ -a, Dt = Dt - & 
Hence, with W = DtW and W = DtW, we can write the equation (2.60) as 

(2.62) LiW" = - g^''db{idcP)W'' - q^) - B^iW, W) 
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where qi is chosen so the divergence of the second term on the right vanishes and 

(2.63) gabB\W, W) = -{DtQac - ^ac - 'i^9ac)W^ + {^{Dtgac " C^ac ' ^gac)W' " 9„,go 

Here g'o is chosen as follows so that divLiW^ = d\vW — Aw B = D"^ A\vW + (i\vW a . But I)j divVF = 
Dl div + 2a bt div + (ct — a^) Ai\W so we must have divi? = 2a Dt divVK — ct^ divVF. Hence go is 
chosen so this is fulfilled and (2.49) follows by writing go = ^2 + 93- (2-54) follows from (2.34) or (2.49) 
It follows from (2.49) that we write Li in the two alternative forms: 

(2.64) gabLlW'' = Dt {gabW'> - {u^ab + ^gab)W'') - da{{d,p)W' - qi) + {DtUJab + ^gab)W' + daqo 

(2.65) gabLiW' = Dt{gabW') - da{id,p)W' - qi) - {coab + agab){W' - dW') - dDtgabW + daqo 

(2.55) and (2.56) follows from these. Finally, we also want to express Lq = $'(.t) is these coordinates. 
In order to to this we must transform the term 5x^dk^^ in (2.33) to the Lagrangian frame. If = 
then {Sx''dk^^)dy°- /dx'^ = l^'^Vc^", where Vc is covariant differentiation, see e.g. [CL], and 
(2.57) follows. □ 

3. The projection onto divergence free vector fields and the normal operator. 
Let us now also define the projection P onto divergence free vector fields by 

(3.1) PC/« = - Apu=dwU, Pu\g^ = 
(Here Aq = K~^da{Kg"'^dbq) . ) P is the orthogonal projection in the inner product 

(3.2) {U,W) = [ gabU'^W'^Kdy 

Jn 

and its operator norm is one: 

(3.3) \\PW\\ < \\W\\, where \\W\\ = {W,W)^/^. 

For a function / that vanishes on the boundary define AfW" = g'^^Aj^Wb, where 

(3.4) AfWa = -da{{dJ)W'^-q), A{{dJ)W'^-q)=0, q\Q^ = 0, 

i.e. AfW is the projection of —g°'''db{{dcf)W'^). This is defined for general vector fields but it is only 
symmetric in the divergence free class. We have 

(3.5) {U,AfW)=f naU''{-dJ)WdS, if divU = divW = 0, 

where n is the unit conormal. f\g^ = then ~dcf\g^ = {—^Nf)nc- It follows that ^4/ is a symmetric 
operator on divergence free vector fields, and in particular, the normal operator in (2.50) 

(3.6) A = Ap 

is positive since we assumed that — VjvP > c > on the boundary. We have 

(3.7) \{U,AfW)\<\\VNf/VNP\\L^{dn){U,AU)^'^{W,AW)^'\ if divt/ = divl^ = 0. 
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Since the projection has norm one it follows from (3.4) that 

(3.8) \\AfW\\ < \\d^f\\L^in)\\W\\ + mU^^n^WdWl 

Note also that Af acting on divergence free vector fields by (3.5) depends only on '^n/Iqq, i-e. = Af 

if Vjv/Iq^ = ^nIIqq- We can therefore replace / by the Taylor expansion of order one in the distance 
to the boundary in polar coordinates multiplied by a smooth function that is one close to the boundary 
and vanishes close to the origin. It follows that 

(3.9) ||^/W^||<<:7 5^||Vjv5/|Uoc(an)||W^||+C||Viv/||Loo(9n)(||5Ty|| + ||W^||), if divT^=0, 

ses 

where iS is a set of vector fields that span the tangent space of d^, see section 4. 

In order to prove existence for the linearized equations we in [LI] replaced the normal operator A 
by a smoothed out bounded operator that still has the same positive properties as A and commutators 
with Lie derivatives, and which also has vanishing divergence and curl away from the boundary. This 
makes possible to pass to the limit and obtain existence for the linearized equations. The smoothed out 
normal operator is defined as follows. Let p = p(d) be a smooth out version of the distance function 
to the boundary d{y) = dist(y, 50) = 1 — |?/| in the standard Euclidean metric Sijdy^dy^ in the y 
coordinates, > 0, p{d) = d, when d < 1/4 and p{d) = 1/2 when d> 3/4. Then we can alternatively 
express Af as 

(3.10) AfWa = -da{{f/p){dep)W'-q), A{{f/p)id,p)W''-q)=0, g|a^ = 

Let x{p) be a smooth function such that x' ^ 0, xip) = when p < 1/4, x{p) = 1 when p > 3/4. Af is 
unbounded so we now define an approximation that is a bounded operator: A'yW"' = g°'''A^fWi,, where 

(3.11) AjWa = -Xeda{{f/pmp)W'^)+daq, Aq = da{g''^ KXedl,{{f / p){d,p)W'^)) , g|^^=0 

where Xe{p) = x{p/^)- We have 

(3.12) {U,A}W)= J^{f/p)x'eidaP)U''{dcp)W^Kdy, if divt/ = divW^ = 0, 

from which it follows that A'j is also symmetric. And in particular A^ = A^ is positive since we assumed 
that p > 0, at least close to the boundary. We have 

(3.13) \{U,A}W)\<\\f/p\\Lo.^n\n^,,){U,A'U)^/^{W,A'Wy/^, if divU = divW = 0, 
where = {y € ^l;d{y,dQ) < e}. It also follows from (3.12) that another expression for A^f is 

(3.14) A}Wa = {f/p)x'e{daP) {dep)W'^-daq, Aq = K-'da{Kg''\f/p)x'e{dbP) idcP)W^) , q\g^= 
acting on divergence free vector fields. Furthermore, by (3.12) 

k 

(3.15) \\D'^A'W\\r<Cs^\\DiW\\r, where \\W\\r = ^ \\d^W{t,-)\\L-(n) 

j=0 \oi\<r 
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Let us also define the projected multiplication operators Mg with a two form (3 by 

(3.16) M^Wa = EWabW'') 

Since the projection has norm one it follows that 

(3-17) \\M^W\\ < ||/3||oo||W^|| 

Furthermore we define the operator taking vector fields to one forms 



(3.18) GWa = M^Wa = PigabW) 

Then G acting on divergence free vector fields is just the identity /. 

Let Li be the modified hnearized operator in (2.49) and let W = DtW = DfW + (divV)W = 
k~^Di;{kW), W = DfW. We want to prove existence of a solution W to 

(3.19) LiW = W + AW - BqW - BiW = F, w\^^^^ = W\^^^ = 

for general vector fields F that are not necessarily divergence free. To do this we first subtract of a 
vector field Wi that picks up the divergence and then solve (3.19) in the divergence free class. Let us 
decompose a vector field into a divergence free part and a gradient using the orthogonal projection: 

(3.20) W = Wo + Wu Wo = PW, W^f=5"'56gi, Qi\an = ^- 

Then if gab = DtQab, where Dt = Dt - ct, we have daDtqi = Dt{gabW^) = gabW^ + gabW^ and 
daD^qi = gabW^ + 2gabW^ + gabW^, where gab = D^gab- Hence 

(3.21) = g'^'dbD^q, - 2g-''gbcW^ - g'^'ghcW^ 

Since D^qi\g^ = and the projection of a gradient of a function that vanishes on the boundary vanishes 

(3.22) PW^ = B2{Wi,W,r, where B2{Wu Wi)" = - P {2g"'' gbcW^ + g'^^gbcW^) 
Since divl^/Q = it follows that divWo = divH^o = and hence by Lemma 2.5 

(3.23) PLiWo = LiWq = Wq + AWq - BiWq - BqWo 

(3.24) PLiWi = AWi - BiiWi - BqiWi 
where 

(3.25) SiiVF"= P5iVF"+2P(rt6cW^'=) BoiW''= PBoW''+P{g''^gbcW). 
Hence projection of (3.19) gives 

(3.26) LiWo = -PLiWi + PF = -AWi + B^Wi + BoiWi + PF, 
Here, by (2.54) 

(3.27) W^=g''%qi, Aqi = ip, g^f, = 0, 
where 

(3.28) D'^ip + aip = divF. 
By (3.23)-(3.24) we also have 

(3.29) {I-P)LiWo = 

(3.30) (/ - P)LiWi = Wi- B2{Wi, Wi) -{I- P)BoW + (/ - P)BiWi 
Summing up, we have proven: 
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Lemma 3.1. Suppose that W satisfies L^W = F. Let Wq = PW , Wi = {I - P)W , Fq = PF and 
Fi = {I- P)F. Then 



(3.31) LiWo = -Fo - AWi + BnWi + Bq^Wi 

(3.32) =F^+ B2{Wr, W^) + {I - P)BoWi + {I - P)BiWi 

where Bqi and Bn are given by (3.25), B2 is given by (3.22) and Bq, Bi is as in (2.51), (2.52). 
Furthermore 

(3.33) D'^ div Wi + 'a div Wi = div F 



We now find a solution of (3.19) by first solving tlie ordinary differential equation (3.28) and then 
solving the Dirichlet problem for qi and defining Wi by (3.27). Finally we solve (3.26) for Wq within 
the divergence free class. This gives existence of solutions for (3.19) for general vector fields F once we 
can solve it for divergence free vector fields. However, we also need estimates for (3.19) that do not 
loose regularity going from F to in order to show existence also for the linearized equations (2.57): 

(3.34) LoW = L{W - BzW = F, = = 0, 

by iteration. It seems like there is a loss of regularity in the term —AWi in (3.26). However, curl^l^^i = 
and there is an improved estimate for (3.19) when divF = and curlF = 0, obtained by differentiating 
with respect to time and using that an estimate for two time derivatives also gives an estimate for the 
operator A through the equation (3.19). We can estimate any first order derivative of a vector field in 
terms of the curl, the divergence and the normal operator A and there is an identity for the curl. 

Let us now also derive the basic energy estimate which will be used to prove existence and estimates 
for (3.19) within the divergence free class: 

(3.35) W + AW = H, wl^^ = W\^^^ = 0, dwH = 

where A is the normal operator or the smoothed version. For any symmetric operator B we have 

(3.36) — (VF, BW) = — KWBWa dy = 2{W, BW) + {W, BW) 
at at Jq 

where W = K~^Dt{KW) and B is the time derivative of the operator B considered as an operator from 
the divergence free vector fields to the one forms corresponding to divergence free vector fields: 

(3.37) BW = P{g''\DtBWb - BWb)) , BWb = QbcBW^, 

see section 4. The projection comes up here since we take the inner product with a divergence free 
vector field in (3.37). Let the lowest order energy Eq = E{W) be defined by 

(3.38) E{W) = {W, W) + {W, {A + I)W) 

Since {W,W) = {W,GW), where G is the projection onto divergence free vector fields given by (3.18), 
it follows that 



(3.39) 



^0 = "^{W, W + {A + I)W) + {W, GW) + {W, (A + G)W) 
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In particular it follows from (3.4) or (3.10) respectively (3.16) and (3.18) that 



(3.40) Af = Ai, G = Mg, where f = nDtiK'^ f) and g = KDt{K-^g). 

In fact the time derivate of an operator, as defined by (3.37), commutes with the projection since 
DtdaQ = daDtq, where Dtq\g^= if Q\gQ= 0, and the projection of the gradient of function that 
vanishes on the boundary vanishes. It therefore follows from (3.7) or (3.12) and (3.17) that 

(3.41) \{W,Aw)\<\\p/p\\^{W,AW), \{W,GW)\<\\g\\^{W,W) 

The last two terms in (3.38) are hence bounded by a constant times the energy so it follows that 

(3.42) |i;o| < /^(2||iJ||+c/^), c= ||p/p|U + ||5lloo + 2 

from which a bound for the lowest order energy follows. 

Similarly, we get higher order energy estimates for vector fields that are tangential at the boundary, 
see section 10. Once we have these estimates we use that any derivative of a vector field can be 
bounded by tangential derivatives and derivatives of the divergence and the curl, see section 5. The 
divergence vanishes and we can get estimates for the curl as follows. Let Wa = QabW^, Wa = gabW^ 
and iba = gabW^ ■ Then DtWa = gabW' + Wa and DtWa = gabW^ + tUa where gab = Dtgab = nDtingab)- 
Since 

(3.43) iv + AW = H, H = BoW + BiW + F 
where curl^PF = it follows that 

(3.44) |Acurlu;| + \Dtcmlw\ < C{\dW\ + \W\ + \dW\ + \W\ + |curlF|) 

Note that the estimate for the curl is actually very strong. The higher order operator A vanishes so 
there is no loss of regularity anymore and furthermore the estimate is point wise. This crude estimate 
suffices for the most part. However, there is an additional cancellation, whereas one would not need to 
assume estimate for \dW\ in the right hand side of (3.41). The improved estimate is for Wa replaced by 
Wa = Wa — ^abW', where Uab = daVb — ObVa- It follows from Lemma 2.5 that 

(3.45) \Dtcmlw\ + \Dtcmlw\ < C(| curlw}| + |ai^| + |I^| + | curlF|), \cuA{w-w)\ < C{\W\ + \dW\) 

4. The tangential vector fields. Lie derivatives and commutators. 

Following [LI], we now construct the tangential vector fields, that are time independent expressed in 
the Lagrangian coordinates, i.e. that commute with Df. This means that in the Lagrangian coordinates 
they are of the form {y)d / dy"- . Furthermore, they will satisfy, 

(4.1) daS'' = 0, 

Since J7 is the unit ball in R" the vector fields can be explicitly given. The vector fields 



(4.2) 



y^d/dy'' - y'd/dy'^ 
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corresponding to rotations, span the tangent space of the boundary and are divergence free in the 
interior. Furthermore they span the tangent space of the level sets of the distance function from the 
boundary in the Lagrangian coordinates 

(4.3) d{y) = dist{y,dn) = l-\y\ 

away from the origin y 7^ 0. We will denote this set of vector fields by 5o We also construct a set 
of divergence free vector fields that span the full tangent space at distance d{y) > do and that are 
compactly supported in the interior at a fixed distance do/ 2 from the boundary. The basic one is 

(4.4) h{y\ yn[fiy')9'iy')d/dy' - r{y')g{y^)d/dy^) , 

which satisfies (4.1). Furthermore we can choose f,g,h such that it is equal to d/dy^ when |?/*| < 1/4, 
for i = I, ...,n and so that it is when \y^\ > 1/2 for some i. In fact let / and g be smooth functions 
such that f(s) = 1 when |s| < 1/4 and f(s) = when |s| > 1/2 and g'{s) = 1 when |s| < 1/4 and 
g{s) = when \s\ > 1/2. Finally let h{y^, y") = f{y^) ■ ■ ■ f{y"')- By scaling, translation and rotation 
of these vector fields we can obviously construct a finite set of vector fields that span the tangent space 
when d> do and are compactly supported in the set where d> do/2. We will denote this set of vector 
fields by Si. Let S = So U Si denote the family of tangential space vector fields and let T = S U {Dt} 
denote the family of space time tangential vector fields. 

Let the radial vector field be 

(4.5) R = y^d/dy". 
Now, 

(4.6) = n 

is not but for our purposes it suffices that it is constant. Let TZ = SU {R}. Note that TZ span the full 
tangent space of the space everywhere. Let U = S U {R} U {-D*} denote the family of all vector fields. 
Note also that the radial vector field commutes with the rotations; 

(4.7) [R, S] = 0, SeSo 

Furthermore, the commutators of two vector fields in So is just it another vector field in So- Therefore, 

for i = 0, 1, let 7^i = 5^ U {R}, Ti = SiU {Dt} and Ui = SiU {R} U {Dt}. 

Let us now introduce the Lie derivative of the vector field W with respect to the vector field T; 

(4.8) CtW = TW - {dcT'')W'' 

We will only deal with Lie derivatives with respect to the vector fields T constructed above. For those 
vector fields T we have 

(4.9) [Dt,T], and [A,>Cr]=0 
The Lie derivative of a one form is defined by 



(4.10) 



Cxaa = TUa + {daT'')ac, 
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The Lie derivative also commute with exterior differentiation, [jCt, d] = so 



(4.11) CrdaQ = daTq 

if g is a function. The Lie derivative of a two form is given by 

(4.12) CrPab = Tpab + (.daT^)Pcb + idbT-)Pac 

Furthermore if w is a one form and cml Wab = dwab = daW^ — d^Wa then since the Lie derivative 
commutes with exterior differentiation: 

(4.13) CT^MxlWab = CUllCTWab 

We wiU also use that the Lie derivative satisfies Leibniz rule, e.g. 

(4.14) Cria^W) = {CTa,)W' + a^CrW, CriPacW) = {CTfiac)W' + (iacCrW^. 
Furthermore, we will also treat Dt as if it was a Lie derivative and set 

(4.15) Cd, = A 

Now of course this is not a space Lie derivative. It can however be interpreted as a space time Lie 
derivative restricted to the space components. What we use is that it satisfies the same properties 
(4.9)- (4. 14) as the other Lie derivatives wc are considering. The reason we want to call it C^t is simply 
a matter of that we will apply products of Lie derivatives and Df and since they behave in exactly the 
same way it is more efficient to have one notation for them. 

The modification of the Lie derivative 

(4.16) CuW = CuW + (divC/)W, 
preserves the divergence free condition: 

(4.17) diY CuW = UdbjW, where Uf = Uf + {diYU)f. 

if / is a function. (4.16) is invariant and (4.17) holds for any vector field U. However, since we are 
considering Lie derivatives only with respect to the vector fields constructed above and only expressed 
in the Lagrangian coordinates it is simpler to use the modification 

(4.18) CuW = K-^CuiKW) = CuW + {Ua)W, where a = In/t 

Due to (4.1), divS" = K~^da{KS"') = Sa, if S is any of the tangential vector fields and divi? = Ra + n, 
if R is the radial vector field. For any of out tangential vector fields it follows that 

(4.19) divCuW = UdwW, where Uf = Uf + {Ua)f = K-^U{Kf). 

This has several advantages. The commutators satisfy [Cu,Ct] = C[u,t]: since this is true for the usual 
Lie derivative. Furthermore, this definition is constant with our previous definition of D^. 
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However, when applied to one forms we want to use the regular definition of the Lie derivative. 
Also, when applied to two forms most of the time we use the regular definition: However, when applied 
to two forms it turns out to be sometimes convenient to use the opposite modification: 

(4.20) CrPah = jC-Tpab - {Ua)pab, 

We will most of the time apply the Lie derivative to products of the form aa = PabW': 

(4.21) CriPabW) = {CT/3ab)W' + pabCrW 

since the usual Lie derivative satisfies Leibniz rule. Using the modified Lie derivative we indicated 
in [L2] how to extend the existence theorem in [LI] to the case when k is no longer constant, i.e. 
Dta = divV ^ 0. This will be carried out in more detail here. 

Let U = be some labeling of our family of vector fields. We will also use multindices 

/ = (zi, ...,ir) of length |/| = r. Let = Ui^ ■ ■ ■ Ui^ and jCfj = jC-m^ ■ ■ ■ ^Ui^, where jCu is the Lie 
derivative. Similarly let f = Ui^ ■ ■ ■ UiJ = K-W\nf) and t\jW = tu,^ ■ ■ ■ ^Ui,W = k-^C{j{kW), 

where Cu is the modified Lie derivative. Sometimes we will also write Cfj, where U E Sq or I E Sq, 
meaning that Uif, G So for all of the indices in /. 

We will now calculate commutator between Lie derivatives and the operator defined in section 3, i.e. 
the normal operator and the projected multiplication operators. It is easier to calculate the commutator 
with Lie derivatives of these operators considered as operators with values in the one forms. The one 
form w corresponding to the vector fields W is given by lowering the indices 

(4.22) Wa = W,= gabW^ 

For an operator B on vector fields we denote the corresponding operator with values in the one forms 
by 5. These are related by 

(4.23) BWa = 9abBW\ BW = g'^'B, 

Most operators that we consider will map onto the divergence free vector fields so we will project the 
result afterwards to stay in this class. Furthermore, in order to preserve the divergence free condition 
we will use the modified Lie derivative. If the modified Lie derivative is applied to a divergence free 
vector field then the result is divergence free so projecting after commuting does not change the result. 
As pointed out above, for our operators it is easier to commute Lie derivatives with the corresponding 
operators from the divergence free vector fields to the one forms. Let Bt be defined by 

(4.24) BtW" = P{g'''' {CrRWh - Bf^CrW)) 

In particular if i? is a multiplication operator B_^W = P{PabW') = PabW' — daQ, where q vanishes on 
the boundary is chosen so that div BW = then 

(4.25) jOtB^W = PabJC-rW' + {jC-TPab)W'' + daTq 

and if we project to the divergence free vector fields then the term daTq vanishes since if T is a tangential 
vector field then Tq = as well. It therefore follows that Bt is another multiplication operator: 



(4.26) 



BrWa = P{{£TPab)W') 
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In particular, we will denote the time derivative of an operator hy B = B^t and for a multiplication 
operator this is 

(4.27) BW = Bd,W = P{{DtPab)W^) 

If B maps on to the divergence free vector fields 

(4.28) ItBW" = It{9''''B^W) = {jC,Tg"''')B^W + q^^CtB^W 

Here Ctq"'^ = —g°''^g^'^^T9cd- If B maps onto the divergence free vector fields then CtB is also 
divergence free so the left hand side is unchanged if we project: 

(4.29) CtBW = -P{g''\CTgbc)BW^) + Pig^'^CxB^W - B^CtW)) + BCtW 

By (4.26) applied the G„5 = PigabW'') we see that GtW = P{{g''^ CTgbc)W'')) so the first term in the 
right of (4.29) is GtBW^. The second term is by definition (4.24) BtW so we get 

(4.30) CtBW = BCtW + BtW - GtBW 

The most important property of the projection is that it almost commutes with Lie derivatives with 
respect to tangential vector fields. If Pua = Ua — daPu then 

(4.31) PC-TPUa = PC-TUa 

since CTdaPu = daTpu vanishes when we project again since Tpu vanishes on the boundary. We have 
just used this fact above. We have already calculated commutators between Lie derivatives and the 
multiplication operators so let us now also calculate the commutator between the Lie derivative with 
respect to tangential vector fields and the normal operator. Recall that the normal operator is defined 
by AfW°- = g''''AfWb, where 

(4.32) AfWa = -d^{idJ)W'=-q), A{{dJ)W--q)=0, q\9^ = 

and / was function that vanished on the boundary. Since the Lie derivative commutes with exterior 
differentiation it follows that 

(4.33) jC-TAfWa = -daiidj)^^ + {d.TDW + {d.Ta)!^ - Tq) 

Since q vanishes on the boundary it follows that Tq also vanish on the boundary and so does {dcTa)fW'^. 
Therefore the last two terms vanish when we project again so we get 

(4.34) P(5«^£T^^W^b) = P{g''''A^tTWb) + P [g"'' At f^b) 

Let us now change notation so yl = ^p, where p is the pressure. Then we have just calculated At 
defined by (4.24) to be At = Af^, i.e. 

(4.35) AT = Afp, if A = Ap 

In particular, if T = is the time derivative we will use the notation A = Ad^ which then is 

(4.36) Aw = Ad,W = A^^pW 
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Exactly the same formulas hold for ^y-. By (3.14) 

(4.37) A}Wa = {f/p)Xe{daP) {d,p)W'^-daq, Aq = K-^da{Kg'^''{f/p)XeidbP) (5cP)W^), q\9^=0 

where p = p{d), d{y) = dist (y, 90). It follows that Tp =, if T e %. Furthermore S e Si = S \ So 
vanishes close to the boundary when d{y) < do/2 and Xe = when d{y) > £ so it follows that 

(4.38) CrAjWa = {{ff)/p)xUdaP) {dcP)W^- {f / p)x'e{daP) {P^pYCtW^- d^Tq. 
Hence 

(4.39) P{g''''CTA)Wk) = P {g'^'^ A) CtW^) + P{g^''MrfWb) 

We can now also calculate higher order commutators: 

Definition 4-1- If T is a vector fields let Bt be defined by (4.24). If T and S arc two tangential vector 
fields we define Bts = {Bs)t to be the operator obtained by first using (4.24) to define Bs and then 
define {Bs)t to be the operator obtained from (4.24) with Bs in place of B. Similarly if = 5*^ • • •S'*'' 
is a product of r = |/| vector fields then we define 

(4.40) Bj = {-..{Bsn)---)s.. 

If is a projected multiplication operator BW"" = P^g""^ l^bcW'^^ then 

(4.41) BiW = P{g''\Cir(5bc)W^) ■ 
In particular if GW = P{g''^gbcW^) then 

(4.42) GiW = P{g''\£'r9bc)W'^). 
If A is the normal operator then 

(4.43) AiW = P{g''^db{{dcf^p)W'')) 

With Bt as in (4.4) we have proven that if B maps onto the divergence free vector fields then 

(4.44) CtBW = BWt + BtW - GtBW, Wt = CtW 
Repeating this, gives for a product of modified Lie derivatives: 

(4.45) C^rBW = c{'-^''Gi, ■ ■ ■ Gi.Bi^Wi, Wj = C^W 

where the sum is over all combinations of / = /i + ... + 7^, and c/^ "^'' are some constants, such that 
^7i...7fc = 1 if /-^ -)- = I. Let us then also introduce the notation 

(4.46) G/^^^ = c/-^'=Gi3 

where the sum is over all combination such that I3 + ...1^ = I — Ii — I^- With this notation we can 
write (4.41) 

(4.47) C^tBW = Gj'^'Bi^Wi^ 
where again G/'-^^ = 1 if Ji + /2 = /. Also let 

(4.48) G/^-^'==0, if I2 = I, and Qh-ik ^ Qh-h ^ otherwise. 
Then we also have 

(4.49) C^tBW = BWi + Gj'^^Bi^Wi^ 
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5. Estimating derivatives of a vector field in terms of the curl, 
the divergence and tangential derivatives or the normal operator 

The first part of the lemma below says that one can get a point wise estimate of any first order 
derivative of a vector field by the curl, the divergence and derivatives that are tangential at the boundary. 
The second part say that one can get estimates with a normal derivative instead of tangential 
derivatives. The last part says that we can get the estimate for the normal derivative from the normal 
operator. The lemma is formulated in the Eulerian frame, i.e. in terms the Euclidean coordinates. 
Later we will reformulate it in the Lagrangian frame and get similar estimates for higher derivatives. 

Lemma 5.1. Let Sf be a vector field that is equal to the normal M at the boundary dVf and satisfies 
\M\<1 and\dM\<K. Let q'^ = 6'^ -M'N^. Then 

(5.1) < C{q^'5'^dkPi dipj + I curl/3|2 + | div/3p) 

(5.2) / \d(3fdx <C [ {5^^M^M^dil3k djf3i + | curl/?|2+ | div/3p+ K^\p\'^) dx 

Suppose that d^^aj is another vector field that is normal at the boundary and let A(5i = di{akf5^ — q) 
and q is chosen so that divA/3 = and g|an = 0. Then 

(5.3) f 5'^akaidiP''djpUx<C [ {5'^APiApj + \af{\curlpf+\divPf)+\daf\pf)dx 

Proof. (5.1) follows from the point wise estimate 

(5.4) 6'^S''^WkiWij < C{S'^q'''wkiWij + \w\^ + (tiwf) 

(5.5) S'^S'^'wkiWij < C{N'N^5^'wkiWij + {q'^q'^' - q"'q^')wkiWij + \wf + (tr wf) 

where Wij = Wij — wji is the antisymmetric part and triy = S'^^Wij is the trace. To prove (5.4)-(5.5) 
we may assume that w is symmetric and tracclcss. Writing (5*' = q^^ + M'^M^ we see that (5.4) for 
such tensors follows from the estimate Af^Af^ Af^Af''WkiWij = {N'^N^WkiY = {(f^WkiY < nq'^-' q^^WkiWij- 
(This says that {tT{QW)f <niT{QWQW) which is obvious if one writes it out and use the symmetry.) 

(5.5) follows since {5^^ q^^ - M'-Sf^ 6^^)wkiWij = {q"-^ q^^ - M""]^^ N'^N^)wkiWij = {q^^'^^ - q''''q^^)wkiWij. 
(5.2) follows from (5.5) and integration by parts using that the boundary terms vanish, since we assumed 
that M = M there, and that {q^^q''^ - q'^q^^)^i dkdjpi = 0: 

(5.6) / {q'^q'^' - q'\=')dkPi dj/3i dx = - [ dkiq^'q"' - q'\'') f5i d^f5i dx 

JVt JVt 

We have APi = {diak)fi^ + akdifi^ — diq so to prove (5.3) we must estimate ||9(7||7;,2. Since = 
diAfi^ = A{akP^) - Aq it follows that Aq = A{ak(3^) = 25^ ( (5* )/?'=) + afcA/?'^ - {Aak)P^ and 
akA(3^ = ai(aMiv/3 + afeCurl/3*'=) - diva div/3 - {dkai)d''P' + {dkai)d'P'', and hence A(afc/3'=) = 

di (2{d'ak)P'' + a' div/3 + Uk cmip''' - diva - curla\ /J'^) . It follows that 

(5.7) / \dqfdx = -[ qAqdx = - f qdi{a' dW p + akCuiip''' + {d'ak + dka')p'' - diva P') dx 
Jci Jn Jn 

23 



and integrating by parts again gives ||9g||i2 < C(|| |a|div/?||L2 + || |a| curl/3||i2 + || |5a|/3||i2). □ 

Definition 5.1. For V any of the family of vector fields introduced in [LI] and for (3 a two form, a one 
form, a function or a vector field we define 

(5.8) \p\]!= Yl \^'uP\, [P]^= E 

|/|<r,7eV ri + ...rk<r,ri>l 

and [/3]o = 1. Furthermore let 

(5-9) IPlr = E 

\oi\<r 

If /3 is a function then CuP = U {3 and in general it is equal to this plus terms proportional to (3. 
Hence (5.8) is equivalent to just the sum X^|7|<^ (3\ . In particular if TZ denotes the family of 

space vector fields then is equivalent to \P\r with a constant of equivalence independent of the 
metric. Note also that if (3 is the one form = daQ then jC,fjj3 = dU^q so \dq\^ = Yl\i\<r lev \dU^q\. 

Definition 5.2. Let ci be a constant such that 

n 

\dxldyf + \dyldxf < cl E i\9ab\ + Iff"'!) < ncj, 

a, 6=1 

and let Ki denote a continuous function of ci. 

We note that the second condition in (5.10) follows from the first and the first follows from the 
second with a larger constant. We remark that this condition is fulfilled initially since we are composing 
with a diffeomorphism. Furthermore, for solution of Euler's equations, divF = 0, so the volume form k 
is preserved and hence an upper bound for the metric also implies a lower bounded for the eigenvalues 
and an upper bound for the inverse of the metric follows. 

In what follows it will be convenient to consider the norms of CfjW = k,~^C(j{k,W) if is a vector 
field and of £.(jg = KjCfj{K^^g), if g is the metric. The reason for this is simply that div {jCfjW) = 

dwW and CljCmlw = curl{Cfjw) and when we lower indices Wa = gabW^ = {K~^gab){KW') and 
apply the Lie derivative to the product we get JCuWa = {jC.ugab)W' + gabtuW^. 

Lemma 5.2. Let W he a vector field and let Wa = gabW' be the corresponding one form. Let k = 
del {dx / dy) = ^/detg . Then 

(5.11) \k\ + \k-^\<Ki, \U^K\ + \U^K-^\<Kic^'-^''\U^'g\---\U^''g\ 

where the sum is over all It + ...Ik = I. 

With notation as in Definition 5. 1 and section 4 we have 

r-l 

(5.12) \KW\f < Ki{\cuvlw\f_i + |KdivW|^_i + \KW\f + J2\9/'^\'^-s\'^W\f) 

s=0 
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(5.10) 



We also have 

r 

(5.13) \kW\'^ < ^1 E (|curlii;|^_i_, + \KdiyW\'^_,_, + 

where for s = r we use the convention that | curli(;|};^j = |«;divVF|Yi = 0. Furthermore (5.12)-(5.13) 
holds without the factors n and l/n, i.e. 

r 

(5.14) \W\f <K^Y. [5]?(|curlu;|?_i_, + \db,W\f_^_, + \W\f_,), 
(5.12)-(5.13) also holds for the vector field W replaced by a one form w, i.e. 

r 

(5.15) \w\f <K^Y, [5]?(|curlu;|^_i_, + \d\^W\f_,_, + \w\f_,), 
Moreover, the inequalities (5.12)-(5.15) also hold with {TZ,S) replaced by {U,T). 

Proof Ua = luK = (lndet£f)/2 then Ua = tYCug/2 = g^^CuQab/'^ and Cug"^ = -g'^'g^'^CuQcd- An 
easy consequence of Lemma 5.1, see [LI], is: In the Lagrangian frame we have, with Wa = W ^ = gabW', 

(5.16) \CuW\<Ki(\cnrlW\ + \diyW\ + ^s^s\^sW\ + [g]i\W\y UeU, 

(5.17) \CuW\<Ki(^\cmlW\ + \diYW\ + J2Ter\^TW\ + [g]i\W\y U eU. 

where [^Ji = 1 + l^^l. Furthermore 

(5.18) \dW\ < K, {\jC,rW\ + Eses\^sW\ + [g]i\W\^ 

When d{y) < do we may replace the sums over S by the sums over Sq and the sum over T by the sum 
over Tq. In [LI] this was proven for Cu replaced by jCu, but the difference is just a lower order term. 
We claim that 

(5.19) Yl \^hW\<K^ Yl (Icurl^^I^I + jdiv^^I^j + [g]^\jC,ijW\) + ^ \^sW\ 
\i\=r,uen \j\=r-i,uen \i\=r,ses 

First we note that there is noting to prove if d{y) > do since then 5 span the full tangent space. 
Therefore, it suffices to prove (5.19) when d{y) < do and with S replaced by Sq and TZ replaced by TZq. 
Then (5.19) follows from (5.16) if r = 1 and assuming that its true for r replaced by r— 1 we will prove 
that it holds for r. If we apply (5.16) to JO-ifW, where |J| = r— 1, we get 

(5.20) \CuC{jW\ < Ki (I curl^^l + | divi:^I^| + ^ \CsCijW\ + [<7]i |i:^I^|) • 

ses 

If C'Ij consist of all tangential derivatives then it follows that \CijC^W\ is bounded by the right hand 
side of (5.19). If £^ does not consist of only tangential derivatives then, since [Cr,Cs] = ^[r,s] = 0) if 
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S G So, wc can write Cs'C.'ijW = C^Cs'W, for some S' G Sq. If we now apply (5.19) with r replaced 
by r — 1 to Cs'Wj (5.19) follows also for r. 

In the Lemma we have Cfj cuvlw = cuvlC^jW which however is different from curl CfjW, We have: 



(5.21) C{jWa = C{j{gabW') = -gabCijW'' + c'j^jy^lCi/W'', where 9ib = ^u9> 



ah 

where the sum is over all Ji + J2 = J and c j^j^ = 1 for | J2I < \J\ cj^j^ = if J2 = J. It follows that 

(5.22) |curl^^-curl/:^u;| <2c\j^{\dg^'\\C^jfW\ + \g^'\\dClfW\), IJ2I < |J|, 

where the partial derivative can be estimated by Lie derivatives. Furthermore, in the Lemma we have 
\U\KdivW)\ = K-^\U^d:rfW\ = k"^! div>C^VF|. (5.13) follows by induction from (5.12). □ 

Definition 5.3. For V any of the family of vector fields introduced in [LI] let 

(5.23) \\w\\^= Yl w^hm, 11^/^11^00= E w^hwwoo 

|/|<r-,7ev |-f|<r-,/ev 

and let 

(5-24) IIW^II, = J2 W^y^W \\W\\r,oo = II^^^IU 

|a|<r |a|<r 

where \\W\\ = ||W||L2(n), \\W\\^ = ||W||L<-(n)- 

It follows from the discussion after Definition 5.1 and (5.11) that ||M^||r is equivalent to ||W||^ 
with a constant of equivalence independent of the metric. As with the point wise estimates it will 
sometimes be convenient to instead use = (kW)||. This in particular true for the 

family of space tangential vector fields S. However instead of introducing a special notation we then 
write . Since k is bounded from above and below by a constant Ki this is equivalent with 

a constant of equivalence Ki. Furthermore, by interpolation < Xi(||^||r + ||W||f) and 

< i^idl^llrllVFll + IIkWII^ ), and our inequalities anyway contain lower order terms of this form, 
so the inequalities below are true either with or without k. 

Lemma 5.3. We have with a constant Ki as in Definition 5.1: 

r-l 

(5.25) IIW^II^ < Ki(||curlw||^_i + ||KdivTy||^_i + + ^ ||5||^_^,oo 

s=0 

and, with the convention that || curlif;||_i + ||divVF||_i = 0, 

r 

(5.26) IIH^II^ <KiY ||5llr-«,oo(||curlu;||,_i + ||/«div^||,_i + WnWf,) 

Proof. This follows from Lemma 5.2 and the interpolation inequalities below in Lemma 6.2. □ 

We can also bound derivatives of a vector field by the curl, the divergence and the normal operator: 
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Lemma 5.4. Let cq > be a constant such that \ Vnp\ > cq > 0, let K2 and K3 be constants such that 
||VArp||Loo(an) < K2 and ^ses l|ViV'S'p||L«>(aQ) < K3. Then 

(5.27) co\\dW\\ < C{\\AW\\ + i^adl curl«;|| + || dwW\\) + {K3 + b]i)||l^||) 

Proof. We want to express (5.2) and (5.3) in the Lagrangian frame. We also want to pick an extension 
of the normal to the interior. If d{y) be the distance to the boundary in the Lagrangian frame, since is 
the unit ball this is just 1 — \y\. Let Xiid) be a smooth function that is 1 close to and when d > 1/2. 
If Uc = dcd then Uc = Uc/\/ g"-^UaUb is the unit conormal at the boundary and fic = Xi{d)nc defines an 
extension to the interior and N°' = g°'^ni, is an extension of the unit normal to the interior. Similarly, 
by the remarks in section 3, the normal operator only depends on Vjvp restricted to the boundary. Let 
us define = X2{d)fdbd, where / is a function that is equal to N'^dcP = VnP at the boundary and 
extended to be constant along rays through the origin, and X2 is a function that is 1 on the support of 
Xi and when d > 3/4. Then AW = P{g''''db{{dcP)W'')) = Pig^'^dbiacW)) by the remarks in section 
3. Now, in expressing (5.2) and (5.3) in the Lagrangian coordinates partial differentiation becomes 
covariant differentiation so we will pick up a constant coming from the Christoffel symbols, i.e. one 
derivative of the metric [^Ji = 1 + \ dg\. Similarly, one derivative of the normal A''" also gives rise to one 
derivative of the metric. Hence (5.2) and (5.3) become 

(5.28) \\dW\\ < C{\\xiincdW'')\\ + \\cuvlw\\ + \\dwW\\ + b]i||I^||) 
and 

(5.29) \\fX2in,dW'=)\\ < {\\AW\\ + \\fcuilw\\ + ||/divW|| + b]i||/I^|| + || \df\W\\) 

Since |/| > cq and X2 = 1 in a neighborhood of the support of xij the lemma follows. □ 
By Lemma 5.4 we have 

(5.30) co\\dCiW\\ < i^adl cml CjW W + || dwCiW\\ + \\AC-^sW\\ + H-CsV^H) 



where K3 is as in Definition 6.1 and cq as in the physical condition (1.6). Here, the curl of {j0.gW)a = 

9abJ^s^^ t)y (5.22) equal to the curl of CgW plus lower order terms. In particular we see that we can 
get any space tangential derivative in this way so we also get: 

Lemma 5.5. With K3 as in Definition 6.1 we have 

r-l 

(5.31) Co\\W\\r < if3(||curlw;||^_i + ||divW||^_i + + I] ll5l|r-s,oo||'W^IU) 

s=0 

where 

(5-32) \\W\\iA= E 11^4^11 

\i\=s,ies 
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6. Interpolation, the L°° estimates for the pressure 
in terms of the coordinate and the l°° norms. 

Let US now first state the interpolation inequalities that we will use: 

Lemma 6.1. Let 13 he a two form, a function or a vector field. Let ||/9||r he L'^-Sobolev norms and 
\\P\\r,oo is the norms on the unit hall O in R". Then ifO<s<r and j >0 

(6.1) l|/3|U.,oo < C||/3||j;^^/ni/3|i;t,oo 

(6.2) m\s< cmi-'/"- wprJ^ 

For a proof see e.g. [H1,H2], for the L°° norm and [CL], for the L^ norms. ( (6.1) for j > follows 
from (6.1) for j = applied to dy for \a\ < j. )A consequence is: 

Lemma 6.2. With the same assumptions as in Lemma 6.1 we have 

(6.3) ||Q:||j+r-s,oo||^||j+s,oo < (||a|li,oo||/3|li+r,oo + ||/?|li,oo ||a|li+r,oo) 

(6.4) \\l3\\r-s,o.\\W\\s < C{m\o,oo\\W\\r + ||/3||.,oo || W^||o) 



(6.5) ||/l||j+si,oo' ■ ■ ||/fe|li+Sfc,cx) < ll/l ||j,cx) ■ ■ ■ ||/i-l|li,oo||/i||j+si + ...+Sfc,oo||/i+l||j,oo' ' " ||/fc||j,oo 

i=l 

Proof. This follows from using Lemma 6.1 on each factor and the inequality j^^/r iQT--s/r < ^ -|- 5^ e.g. 

(6.6) ii/3ii.-.,ooiiw^ii. < cm\i(Ml:o^^imi-'^'^\\wrj^ 

= C(||/3||o,oo||H^||r)'^'(||/?||r,oo||W^||o)'"'^' < C(||/?||o,oo ||^||. + ||/3||.,oo || W^||o) . 

This proves (6.4). The proof of (6.3) is exactly the same, (6.5) follows from (6.3) by induction. □ 

Let us now introduce some notation to be used in subsequent sections. We will derive tame estimate 
involving the higher norms of the coordinate x with constants that are bounded if some lower norms of 
the coordinate x are bounded: Recall Definition 5.2 of ci: 

n 

(6.7) \dx/dy\^ + \dy/dx\^ < cj, Yl i\9ab\ + Id^'l) < ncj, 

a, 6=1 

and Ki denotes a continuous function of ci . 

Definition 6.1. Let ci be as in Definition 5.2 and for i = 2, 3, 4 let Cj > ci be a constant such that 

(6.8) ||a;||2,oo + ||i||i,oo < C2, 

(6.9) lkll3,oo + ||a;||2,oo + \\x\\l,oo < C3, 

(6.10) lkl|4,oo + ||a;||3,oo + ||^||2,oo < C4 

and let Ki > 1 denote a constant that depends continuously on Cj. 
It now follows from using Lemma 6.2: 
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Lemma 6.3. With Ki as in Definition 5.1 

(6.11) \\dy/dx\\r,oo<Ki\\ 

If di = d/dx'^ = (dy°- 1 dx^^d I dy°- and a = (ai, a„) let 8°' = d^'^ ■ ■ ■ d^" . For any function f we have 

(6.12) \\d"f\\r,oo<M\\f\\r+k,oo + \\x\\r+k,oo\\f\\l,oo), k = \a\. 

Moreover, 



(6.13) 




■{di 


nfn) 


llr.oo < ^l2^ 
i=l 


\\fl\\l,oo^ 


• • l/i-l |l,oo l/i 


|r+l,oo II /i+1 II l,oo ■ 


||/n||l,oo 












+ Ki\\x\ 


1 r+l,oo Wfl II l,oo * 


■'11 fn 1 1 1 , oo ) 




(6.14) 




■{di 


nfn) 


n 

llr.oo < 

i=l 


l/i |l,oo• 


• • l/i-l |l,oo l/i 


|r+2,oo l/i+1 |l,oo' ■ 


||/n||l,oo 












+^^l||x| 


r+2,oo /l l,co' 


II /n II l,oo ) 





Proof. Let A be the matrix dx^/dy"". Using the formula for the derivative of a matrix daA ^ = 
—A~^{daA)A~^ we see that dyA~^ is a sum of terms of the form 

(6.15) A-\d^^A)A-^---{d^^A)A-\ |ai| + ... + jafej = |a| = r 

Since \A~^\ < Cc\ we see from (6.5) with j = that this is bounded by iiri|yl|j.^oo which proves (6.11). 
Now d^d^f is sum of terms of the form 

(6.16) A-\d^^A-') ■ ■ ■ {d^''-^A-')d^''{dyf), + ... + \Pk\ = |7| + |a| = r - 1 + A; 

By (6.5) and what we used proved this is bounded by Ki\\dyf\\r-i+k.<yc^ + Ki\\A\\r-i+k,oo\\dyf\\o.r,c^ 
which proves (6.12). By (6.13) follows from (6.5) with j = and (6.12) with k = 1. (Note by (6.12) 
||9/||o,oo < -?^i||/||i,oo-) Similarly by (6.5) with j = and (6.12) we can bound the left of (6.14) by 



(6-17) ll/ill r+2,oo ~l~ ||^||r+2,oo) ||/2||l,oo ' ' ' ||/n||l,oo 

n 

,cxD + ||2;||2,cxd)||/i||i,oo ■ ■ ■ ll/i— l||l,oo(||/i||r+l,oo + ||^||r+l,oo) 1 1 /i+1 1| l,oo ' ' ' ||/n || l,oo 



i=2 



The first term is of the form in the right of (6.14). The terms in the sum becomes a sum of four terms 
of the form iiri||/ii||2,oo||/i2||r+i,oo multiplied by factors of the form ||/fe||i,oo- Using (6.5) with j = 1 we 
can bound ||/ii||2,oo||^2||r+i,oo < C'||^i||i,oo||^2||r+2,oo + ||^i||r+i,oo||^2||i,oo- This provcs also (6.14). □ 

Lemma 6.4. Let p be the solution of Ap = —{diV^){djV^), where = -Dtx* and let p = Dtp. Then 
for r > 1 we have 

(6.18) ||p||,,oo < i^sdl 

(6.19) ||p||r,oo < K3{\\x\\r,oo + ||a;||r+l,cx. + ||a;||r+2,oo) 
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Proof. We just apply Proposition 7.1 to 

(6.20) = -{diV^)djV\ v' = Dtx\ P|af2 = 

using Lemma 6.3 to estimate the product. (Recall that ||T^||2,oo < -f^s-) Since Dt = dt + V'dk, where 
dt = dt\ ^ , we have 

(6.21) Ap = A{{dt + V''dk)p) = DtAp + iAV'')dkP + 25'^ idiV'')djdkP 
and 

(6.22) A = -{dt + v''dk){{diV^)idjV')) = -2{diV^)idjV') + 2{diV^)idjV'')dkV' 

so 

(6.23) Ap = -2{diVn{djV') + 2{diVn{djV'')dkV' + {AV'')dkP + 2S'^ {diV'')djdkP 

The second part of the lemma now follows from Proposition 7.1 using Lemma 6.3 and the first part of 
the lemma. □ 

Let us now introduce the norms that we will use: 

Definition 6.2. 



(6.24) m,(t) = ||x(t,-)||i+s,oo 

(6.25) ms(i) = \\x{t, ■)\\2+s,oo + \\x{t, ■)\\i+s,oo, 

(6-26) ms{t) = \\x{t, ■)\\3+s,oo + \\x{t, ■)||2+s,oo + Olll+s.oo, 

(6.27) ns{t) = \\x{t, OIU+s.oo + \\x{t, •)||3+s,oo + \\x{t, ■)||2+s,oo, 



We remark that in Definition 5.2 we made an assumption that the inverse of g and dy/dx are 
bounded. This means that ttiq etc are all bounded from below as well. We note that the corresponding 
bounds for the metric Qab = 5ij{dx^ / dy"-){dx^ / dy^) and ujab = {cvLvlv)ab = {dx^ /dy°'){dx^ /dy''){diVj — 
djVi) follows from the bounds for x, x, and x: 

(6.28) Ibllr.oo < Kirrir, \\9\\r,oo + ||<^||r,oo < -f^2?"'r-, llffllr.oo + ||<^||r,oo < -f^S'T^r- 

The proof of this uses the interpolation inequality (6.5) in Lemma 6.2 applied to each term we get when 
we differentiate. In view of the coordinate condition, see Definition 5.1, the same bounds also hold for 
g replaced by the inverse of g. 

Furthermore, we will now prove that the corresponding bounds also for the pressure follows from 
this. We will actually loose a derivative when passing to the bounds for the pressure because we will 
go over Holder spaces, but this does not matter. In Lemma 6.4, we proved that 



(6.29) \\p{t,-)\\r+l,oo<K3mr{t), 

(6.30) \\p{t, •)||r+2,oo + \\p{t, ■)\\r+l,oo < i^3"^,(t), 

(6.31) \\p{t, Ollr+S.oo + \\p{t, •)l|r+2,oo < KsUrit) 
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In particular 



(6.32) |b||2,oo + ||p||l,oo < ^3, \\p\\2,oo < K4 

We will frequently use interpolation, e.g. 

(6.33) rrirrhs < C{mr+s'rho + ?Tio?Tir+s) < K2rhr+sj 

which follows from Lemma 6.1 and the proof of Lemma 6.2 applied to each term we get when multiplying 
any of the expressions (6.24)-(6.27) together. 

We must also ensure that if the physical condition (2.7) and coordinate condition (2.8) holds initially 
they will hold for some small time < t < T, with cq replaced by cq/2 and ci replaced by 2ci. This 
will be proven in section 11, and until then we will just assume that T is so small that4 these conditions 
hold for < t < T. Furthermore, we will also assume that T < cq < 1 since the estimates we will derive 
then will be independent of T and cq . 



7. The L°° estimates for the Dirichlet problem. 
In this section, we give tame Holder estimates for the solution of the Dirichlet problem: 

(7.1) Aq = F, q\g^ = 0. 

Our Holder estimates loose a derivative since we want to use them for integer values. This is not 

important and with an additional loss of regularity, we could have avoided using Holder estimates 
altogether and just gotten the C'^ estimates from the Sobolev estimates, proved in the next section, 
using Sobolev's lemma. Apart from getting estimates for the solution of (7.1) we also need estimates for 
time derivatives and variational derivatives. For this we need to know that the solution of (7.1) depend 
smoothly on parameters if the metric and the inhomogeneous term do. We remark that the coordinate 
condition is critical since it is needed in order to invert the Laplacian. 

One can also use the results in section 5 to get tame estimates for the solution of the Dirichlet 
problem: In fact if we take = g^-^dbq, and Wa = dqQ, then divW^ = Aq and curlu; = 0. Applying 
Lemma 5.2 to W therefore gives: 

(7.2) \W\^<K,J2 [9]f{\^Q\^-i-s + \W\ts), 

s=0 

where for s = r we should interpret |Ag|_i = 0, and 

(7.3) \dq\^ <K,j2 [9]f{\M^-i-s + \dq\f-s) 

s=0 

Therefore it suffices to obtain estimates for tangential derivatives only which is easier because the 
Dirichlet boundary condition is preserved by tangential derivatives. If qlg^ = then S^q\g^ = 0. The 
L°° estimates uses the standard Schauder estimates for the Dirichlet problem. Because we want to have 
the final result in terms of C*^ norms instead of Holder norms these results loose a derivative. 
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Proposition 7.1. If QlgQ = then for r >1 

(7.4) lkl|r,oo < ^3(||Ag||^_i,oo + Ibllr.oollAgllcoo) 

Proof. If we apply Lie derivatives to = g'^^dt,q we get 

(7.5) = g'^'d.S'q + ci^ij'^'^' d,S'^q, g"'' = t'sg-\ Wi = 

and the sum is over all combinations / = Ii + 12, c/^ constants such that c/^ = if I2 = /. Since 
divWj = div>C;^T^ = S^divT^ = S^Aq = K~^S^{KAq) it follows from taking the divergence of (7.5) 
that 

(7.6) A{S'q) = S'Aq - K-'daic^g'^''' d,S'^q), g' = t'^g'^" 

Let ||w||2+a,oo denote Holder norms, see section 17 and Proposition 7.2. By Proposition 7.2 we have 

(7.7) ||5-f5||2+a,oo < i^l(||5^Ag||i,oo +c/^7,(||5^M|l,oo||>5-'^9||2+a,oo + ||^^M|2,oo || ?|| l+a,oo)) 

If we let Mr = l^[/|<r-2 ll'S'^^lb+a.cxj) > 2, = llo'llr+a.oo for r = 0, 1 it follows from Proposition 
7.2 that Mo + Mi < i^3|| Ag||, M2 < K-i\Aq\x,oo and for r > 3 we have: 

r-l 

(7.8) M, <iC3(||Ag||,_i,oo + ^|| 

s=l 

Inductively it follows that 

r-2 

(7.9) Mr < K3[\\Aq\\r-l,oo + ^ ||5f||r-s,oo||Ag||s,oo) < ^3 (|| Ag||r_i,oo + |bl|r||Ag||o,cx)) 

s=0 

where we used interpolation. With I e S, |/| = r — 2 we have hence get from differentiating (7. ) and 
using what we used proved 

(7.10) \\dWi\\ 0,00 < -f'^a (II Ag||r-i,c» + llslkll Ag||o,oo) 

However, once we have bound for the tangential components, the bound for all components in terms of 

these and w Ap. follows from Lemma 5.2. 

Theorem 6.6 in [GT], together with Theorem 8.16, and Theorem 8.33 in [GT] in our setting it reads: 

Proposition 7.2. Suppose that ||^l!>||fe+a,cxD denotes the Holder norms and < a < 1, and k is an 
integer, see section 17. 

(7.11) Ap = g'^'dadhp + K-\daiKg'''))dbP = n-^da {i^g'^'dbp) 
where 

(7.12) ||5"'l|0+a,oo + ||55"'||0+a,oo<A, l^^'l + babj < A 

a,b 
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Suppose that p\gQ = ^- Then 

(7.13) |b||2+a,oo < C'(||p||oo + ||Ap||o+a,oo) 

where C = C{n, a, A, A) and 

(7.14) Iblloo < C||Ap|U 
and ifAp = F + K'^dainG") 

(7.15) Ibll l+a,oo 

< C(||p|U + Ili^lU + ||G|| 

0+a,oo ) 

Furthermore 

(7.16) ||ni;||a,,oo < C'||'"|l7,oo||^^||a,oo, 7>tt, ||'"'y||a,oo < C'(||'"||o,oo||'y||a,oo + ||^^||o,oo||'"||a,oo) 

Note that if we multiply by k then the operator is also in the divergence form that [GT] has in 
Theorem 10.33. Anyway, in our case it is equivalent to a domain in T>t with the the usual metric. 

Let us now prove that the solution of (7.1) depends smoothly on parameters if the metric g and the 
inhomogeneous term F do. Let us assume that the parameter is time t. We have: 

Lemma 7.3. Let ip be the solution of 

(7.17) A(f> = K-^da{K9''''db4>) =F, (t>\g^ = 

where k = \/det g, and g satisfies the coordinate condition (2.8) on [0, T]. Suppose that gab,F S 
C"=([0,T],C°°(n)). Then (f) e C''{[0,T],C'^{n)). 

Proof. Let us write (pt, gt Ft, and = A^^, to indicate the dependence of t. Our initial assumption is 
that guFt& C'=([0, T], C°°(n)). That 

(7.18) At4>t = Ft, 0t|an = O 

has a solution G C°°(n) if -Ft,5t G C°°(ri) and the coordinate condition is fulfilled is well known. 
We win prove that Ft,gte Ci([0,r], C°°(n)) implies that (pt G Ci([0, T], C°°(n)). If this is the case, 
then = Dtcpt satisfies 

(7.19) At^t = Ft-Kt^t, Mdn = ^ 

where At = [A, A] and Ft = DtFt. Since the right hand side of (7.19) is also in C^{[0,T],C°° (Tl)) we 
can repeat argument to conclude that <j)t G C\[Q,T],C°°{Tl)), i.e cpt G C2([0, T], C°°(0)) In general we 
can then use induction to conclude that (pt G C'^{[0,T], C°°{fl)) since 

fc-i 

(7.20) AtD'^cPt = D'^F - ^ cf A^'-^^ Did^t, 

3=0 

where A^^^^ are the repeated commutators defined inductively by A^'^-' = [Dt,A^ "'^-'], A^*^^ = A^. 
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First we will show that Ft,gt G C{[0,T],C°°{n)) implies that (pt G C{[0,T],C°°{n)). We will only 
prove this for i = since the proof in general is just a notational difference from the proof for t = 0. 

(7.21) At{(f>t - (t>o) = Ft - Fo - {At - Ao)^o 

Since the C'^{Cl) or iJ'"(r2) norm of the right hand side tends to as t — for any m and since we 
have uniform bounds for A^^, in Lemma 7.3, it follows that the C™{fl) or H'^ip.) norm of (j)t — 00 
tends to as t ^ for any m. Hence <^t G C{[0,T], C°°{n)). Now, let 4>t be defined by (7.19). By the 
same argument it follows that <j)t G C{[0,T],C°°{il)). It remains to prove that (pt is differentiable. We 
have 

(7.22) At {^t -(t>o- t<po) =Ft-Fo- tFo + {At - Aq - tAo) <^o + * ( A* - Ao)</>o 

Since Ft and gt are differentiable as functions of t it follows that the C"^{Q,) or H'^{^}) norm of the 
right hand side divided by t tends to as t for any m. Since we also have bounds for the inverse 
of A^ that are uniform in t we conclude that any C" or if™" norm of (pt — (po — i^o divided by t also 
tends to as t ^ for any m. It follows that (pt G C^([0, T],C°°{U)). □ 



8. The estimates for the Dirichlet problem. 



In this section, we give tame L^-Sobolev estimates for the solution of the Dirichlet problem: 

(8.1) Aq = F, Q\an = ^- 

We also remark that the coordinate condition is critical in all the estimates in this section since it is 
needed in order to invert the Laplacian A. As pointed out in the beginning of section 7, using the results 
from section 5 it suffices to obtain estimates for tangential derivatives only which is easier because the 
Dirichlet boundary condition is preserved by tangential derivatives. If q^l^j^ = then S^q\g^ = 0. 

Proposition 8.1. Suppose that q is a solution of the Dirichlet problem, qlg^ = 0, and = g"'''df,q. 
Then if r >0 we have 

r-l 

(8.2) \\W\\r<K^^\\g\\r-i-s,oc\\Aq\\,+Ki\\g\\r,oo\\W\\, 

s=0 

and if r > 1 we have 

r-l 

(8.3) \\W\\r + \\q\\r+l < Ki"^ \\g\\r-s,oo\\^q\\s 

s=0 

Furthermore, for i <2 and r > we have 

r-l 

(8.4) \\DiW\\r<K3j2 E \\D^9\\r-s,oo\\DiAq\\s + K3 ^ \\D^g\\r,^\\DiW\\ 

s=0 j+k<i j+k<i 

and if r > 1 we have 



(8.5) 



\DiW\\r + \\Diq\\r+i<K3Y^ E 

s=0 j+k<i 
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Moreover if P is the orthogonal projection onto divergence free vector fields and W is any vector field 
then, for r >0, 



(8.6) \\PW\\r < K,^\\g\\r-s,oo\\W\U 

s=0 

and, for i <2, 

r 

(8.7) \\b\PW\\r<K^Y. E \\Dh\\r-s,oo\\DiW\l 

s=0 j+k<i 

First a useful lemma: 
Lemma 8.2. Suppose that S E S and (l\g^ = 0, and 

(8.8) IsW" = g'^'^dbQ + F". 
Then 

(8.9) ll^s^ll </^i(||divW|| + ||F||) 
Proof. Let Ws = CsW. 

(8.10) / gabWIWlKdy= f W^daQ ndy + f W^gabP'i^dy 

If we integrate by parts in the first integral on the right, using that q vanishes on the boundary we get 

(8.11) [ WSdaqKdy = - [ div {Ws)qKdy 

Jn Jn 

However divVFs = SdivW. Then we can integrate by parts in the angular direction. S = S"'da, 
S = S + dWS so we get J^iSf) Kdy = da{S°-fn) dy = 0, where 5a 5" = 0. Hence we get 

(8.12) [ W^daqndy = [ div {W){Sq)Kdy 

Jn Jn 

Here \Sq\ < Ki\dq\ so it follows that 

(8.13) \\Wsf < /rilldivM^lKllWsll + ||F||) +i^i||W5||||F|| 
and so 

(8.14) llVFsll < -ft^idldivVFll + ||F||) □ 
Proof of Proposition 8.1. If we apply C^g to Wa = gabW^ we get 

(8.15) daS'q = gabW} + c'^''gi,abW}^, Wj = C's^, giab = Cyab 

35 



and the sum is over all combinations I = Ii + I2, c^^^^ are constants such that c/^^^ = if I2 = /. 
If we write = SS"^ , Wj = CsWj, and use Lemma 8.2 we get since divWj = S'''divM^ = S-^Aq = 

(8.16) \\Wi\\ < K^WS-^AqW+K^cl^'^WgiJ^WWiJ 
or if we sum over all of them and use interpolation 

r-l 

(8.17) \\KW\\f <Ki\\KAq\\r-i + KiJ2\\9\\r-s,oo\\KWf,, r>l 

s=0 

We now want to apply Lemma 5.3 to W"' = g^^d^q. Then curliu = and divT^ = Aq so 

r-l 

(8.18) \\W\\r<K^\\nAq\\r-i + K^\\KW\\f + Kr^\\g\\r-s,o^\\W\\s, r>l 

s=0 

We now use (8.17) to replace the term by the terms of the form already in the right hand 

side of (8.18) and we also replace k by 1 which just produces more terms of the same form. Using 
induction in r and interpolation (8.2) follows. 

We also need an estimate for the lowest order term: ||W^||^ = jQg"'^{daQ)idbq) ndy = 
— J^{Aq)q K dy. However, there is a constant depending just on the volume of Q, i.e. k dy, such that 
Ik 1 1 < -f^ill^^ll) see [SY]. Therefore in addition we have 

(8.19) + ll^ll + < iTillAgll, if W'' = g''^dbq. 
Therefore we inductively get from also using interpolation: 

r-l 



^•20) IIkP^II, < i^i^||5||^_,,oo||Ag||,, r>l 



s=0 



We note that we can remove k in the left since doing so just produces lower order terms of the same form. 
This proves the estimate for W in (8.3) and the estimate for q follows from this since = g'^^d^q- In 
fact by (8.15) with replaced by any space vector fields , I & TZ, \\dq\\r < Ki Yl^s=o 11511^-5,00 11^^118 
and by (8.19) we also have an estimate for ||g||. 

It remains to prove the estimate with time derivatives. We can now repeat the argument with i of 
the tangential derivatives being the time derivative, Oj-,^ = Dl and JO.]-)^ = Dl- This gives 

(8.21) daS'Dlq = gabDlW} + c}l'^'^'^{Dl^gi,ab)Drw}^, 

where c/'*i-f2»2 ^ j f^j, (^j^ _^ j^^^^ _|_ ^ ^j^^^ g^^j^ ^j^g^^ jj^j _^ < |7| + z and otherwise. By 
Lemma 8.2 again 

(8.22) \\DiWi\\ < C'||S-^L>^Ag|| +C'c/f^-'^'^||£>^^7,||oopr^/2ll 
where | J| = |/| — 1. Hence 

(8.23) \\DiiKW)\\^ < K^\\Di{i,Aq)\\r-,+K, ||DP(/^-i5)||,_,,<^||L>|(/^T^)||f 
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By Lemma 8.3 below and (8.23) it follows that for z < 2 

(8.24) \\Dl{KW)\\r < Ks(\\Dl{KAq)\\r-, + g)\\r-s,oo\\Di{KW)\U) . 

S'^r, j<i, s+j<r-\-i 

(8.4) now follows from this and interpolation. 

Since DfAq = ADt q + K-'^da{K{btg°-^)dbq) we get 

(8.25) \\ADtq\\ < K2{\\DtAq\\ + \\btg\Woo\\d^ q\\ + || Aff||i,oo||5g||) < K^{\\btAq\\ + \\Aq\\), 

where we used (8.3) with r = 2. Therefore using (8.19) applied to Dtq in place of q we also get an 
estimate for the lowest order norm: 

(8.26) HAWAII + \\dDtq\\ + \\Dtq\\ < K3{\\btAq\\ + \\Aq\\) 

Using this, (8.4) and interpolation gives (8.5) for one time derivative, apart from the estimate 
for By (8.21) with replaced by any space vector fields , I G TZ, \\dDtq\\r < 

Ki{Yll=o IIS'llr-s,oo||^IU + I]s=o IIAS'llr-s.oollW^IU) and by (8.26) we also have an estimate for HA^II- 
Since b^Aq = AD^ q + 2K-^da{K{btg''^)dbDtq) + K-^da{K{b^g''^)dbq) we get 

(8.27) 

\\AD^q\\ < K2{\\b^Aq\\ + ||L>2^||o,oo H^'^H + \\b^g\\i,oo\\dq\\ + || Asllo.oo H^^ A g|| + || Aslli.ooH^A ^zll) 

<iC3(||A'Ag|| + ||AAg|| + ||Ag||) 

where we used (8.5) for i < 1. Therefore we also get an estimate for the lowest order norm: 

(8.28) \\b^W\\ + \\dD^q\\ + \\D^q\\ < K3{\\b^Aq\\ + ||AAg|| + ||Ag||) 

Using this, (8.4) and interpolation gives (8.5) also for two time derivative, apart from the estimate 
for By (8.21) with replaced by any space vector fields , I e TZ, \\dD'^ q\\r < 

^3 I]I=o (Ibllr-s.oollW^IU + II A5llr-s,oo||W^IU + II A^5llr-s,oo || W^||s) and by (8.28) we also have an esti- 
mate for 

It remains to prove the estimates for the projection (8.6)-(8.7). We have W = Wq + Wi, where 
Wo = PW, and Wi = g'^^d^q where Ag = divW^ and g|g^ = 0. Proving (8.6)-(8.7) for r > 1 reduces to 

proving it for r = by using (8.6)-(8.7), since W Aq = div(£|jW) and replacing k by 1 just produces 
more terms of the same form . (8.6)) for r = follows since the projection has norm 1, ||-PVF|| < ||M^||. 
Since the projection of g'^^D^wib = g'^^d^D^q vanishes we obtain from Lemma 8.3 below: 

i-l 

(8.29) ii-pA'^iIi ^^lYl w^r^aw w^i'^iW 

0=0 

Since also PblWo = blWo we have 



(8.30) 



(/ - P)b\Wi = (/ - P)D\W 
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and hence since the projection has norm one 

i-l 

(8.31) \\DiWi\\ + \\DiWo\\ < Ki\\DiW\\ + Ki^\\Di~^ g\\ \\DiWi\\ 

0=0 

Hence for i = 0, 1, 2 it inductively follows that 

i 

(8.32) llA^^oll + llA^^ill < K^^WtWW, i < 2. 

j=o 

Since as before replacing k by 1 just produces more terms of the same form this proves (8.7) also for 
r = 0. (6.4)- (6.5) follows from interpolation. □ 

Lemma 8.3. Let W = g"-^wi, Then 

(8.33) = g'^'Dlw, - J] g'^\br g,,)!)^^ 
Furthermore if W"" = g'^^d^q then 

(8.34) \Dl{KW)\f < K^(\Dl{KdiYW)\f_, + \KW\f + {ng)\f_,\Dl{KW)\f) 

Proof. We have = DKn-^gb^KW) = ^2]^^ Q{Di~^ {K-^gb^))bi{KW''), which proves (8.33). 

(8.34) follows from (5.12) using (8.33) and that fact the curl of Wa = da vanishes to estimate the 
curl of gabbiWK □ 

9. The estimates for the curl. 

We are studying an equation of the general form 

(9.1) W + AW = H, H = B{W,W)+F 

Here is a linear combination of multiplication operators. Here A is the normal operator and it 
projects to the divergence free vector fields even if W is not divergence free. We have curlAVF = and 
div^T^ = so 

(9.2) divW^ = divi7, cutIw = cuilH 

where we defined Wa = gabW^. Now recall that Wa = gabW'' so it follows that 

(9.3) \Dt curlw\ + \DtCUTlw\ < C{\dDtg\\W\ + \dg\dW\ + \dW\ + \cutIw\) 
Similarly 

(9.4) \btdwW\ + lAdivWl < C(|divW^| + |divVt^|) 
Hence 

(9.5) \Dtcm\w\ + |Acurlu;| + |AdivW| + |Adivl^| + \ cutIw\ + |divW| 

< C{\dW\ + \dW\ + \dg\\W\ + \dg\\W\ + |div/f| + |curl^|) 

Since B is of order one and in fact is a multiplication operator it follows that the terms in curl B (TV, W) 
and div B{W, W) are also going to be of the form in the right hand side of (9.5). However, we need to 
take a closer look on what the operator B really is because on the one hand it will give an improved 
estimates and on the other hand we want to find out exactly what the constants above depend on: 
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Lemma 9.1. Suppose that LiW = F, where LiW = W + AW — B{W,W) is given by Lemma 2.4, 
= DtW- and = DjWr Let Wa = gabW\ = 9abW\ Wa = QabWK Then 



(9.6) A CVxlWab = CMTlWab + [gbcW^) " Ob [OacW^) 

(9.7) Dt curlu;„f, = cnrlwab + da [gbcW") - db {gacW"") 

(9.8) curl«;„,, = curlF^^ + cut\B{W, W)ab 

where gab = Dtgab = Dt gab - <^gab and 

(9.9) B,{W, W) = -{gac - iOac - CTgac)W^ + &{gac - u:ac)W^ - daQo 

and LiW = F. On the other hand, if Wa = Wa — {crgab +^ab)W^ and Li is given by (2.54) then 

(9.10) Dt CUllWab = CmlWab + da{{gbc + i^bc + (Tgbc)W) - db{{gac + t^ac + (^gac)W'') 

(9.11) Dtcmlwab = -da{{DtuJbc + ^9bc)W^) + db{{DtuJac + i^gac)W^) +CUrlF„fe 

(9.12) CmlWab = CUllWab + da{i<jgbc + '^bc)W') - db{{<7gac + i^ac)W') 



Proof. The proof uses Lemma 2.5 and the identity A Wa = Dt{gabW'') = gabW^+gabW^ and (2.54). □ 

Now we want to commute with Lie derivatives >C^, since the Lie derivative commutes with the curl: 
CrcvliIw = curli2i^w. 

Using Lemma 5.2 it follows from Lemma 9.1 and Lemma 6.2: 
Lemma 9.2. With notation as in Lemma 9.1 and Definition 6.1 we have 

(9.13) curltullr-i + ll-Dt curltullr-i + || curM||r-i < 2|| curlF||r-i 

r 

+ ^ (||x||^+l_,,oo + \\x\\r+l-s,oo){\\W\\s + \\W\\s) 

We also have 

(9.14) curlu;||r-i + ||-Dt curli(;||r_i < || curlu;||r_i + || curlFHr-i 

T 

+ -f^3 XI (ll^llr+l-s,oo + ||a;||r+l-s,oo + P|| r+l-s,oo) || W^IU 
s=0 

and 

r 

(9.15) I ||curli(;||r-i - || curlii;||r_i | < ^ (||a;||r+i-s,oo + ||a;||r+i-s,oo) II^^IU 

Remark. The difference between on the one hand (9.13) and on the other hand (9.14)-(9.15) is that 
the latter does not require estimates for \\W\\s but instead it requires an extra time derivative of the 
coordinate. However, we do control two time derivatives of the coordinates. 
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10. Existence for the inverse of the modified 
linearized operator in the divergence free class 

We now first want to show that 

(10.1) LiW = W + AW-BoW-BiW = F, Wl^^ = Wo, W\^^^ = Wi, 
has a smooth solution W: 

Theorem 10.1. Suppose that x and p are smooth, p\g^ = and that the coordinate and physical 

condition (2.8) and (2.7) hold for <t <T. Let Li be defined by (2.49) and suppose that Wq, Wi and 
F are smooth and divergence free. Then (10.1) has a smooth solution for <t <T. 

In case, divV = and divF = existence for (10.1) was proven in [LI]. We now want to generalize 
this result to prove existence when divl^ ^ and divF = 0. This is just minor modification of the 
proof in [LI], mostly notational differences, multiplying with k = det (dx/dy) and since divVF = 
ndainW"'). We will just give an outline of the proof. 

First we note that we can reduce to the case with vanishing initial data and F vanishing to all 
orders as f — > 0. In fact, we can get all higher time derivatives by differentiating the equation with an 
inhomogeneous term 

(10.2) b^+^w = Bk{W, ., b^+'^w, dW, db'^W) + b^F, 

where Bk arc some linear functions followed by projections, see (10.8) with / consisting of just time 
derivatives. Let us therefore define functions of space only by 

(10.3) W''^^ = Bk{W°,...,W''^^,dW°,...,dW'')l^^ + b'^Fl^^^, k>Q 
Then 

defines a formal power series solutions at t = 0. What we are doing is just expanding kW in a formal 
power series in t, since Dt{KW) = nbtW. Since divT^'^ = it follows that divT^ = 0. We also note 
that if the initial data are smooth then we can construct a smooth approximate solution W that satisfies 
the equation to all orders as t ^ 0. This is obtained by multiplying the /c*'* term in (10.4) by a smooth 
cutoff xit/^k), to be chosen below, and summing up the infinite series. Here x is smooth x(s) = 1 for 
\s\ < 1/2 and x{s) = for |s| > 1. The sequence £fc > can then be chosen small enough so that the 
series converges in (7"([0, T], iJ*") for any n and m if take (IjW^'^llfc + l)£fc < 1/2. By replacing W by 
W — W and by -F — LiW in (10.1) we reduce to the situation with vanishing initial data and an 
inhomogeneous term F vanishing to all order as t — ^ 0. 

We will therefore assume that initial data in (10.1) vanishes and that F is smooth, divergence free 
and vanishes to all order as t — 0. Then existence of a solution We for the equation where we replace 
the normal operator A by the smoothed out normal operator A^, e > 0, in (10.1) 

(10.5) LlWe = W, + A^W, - BoWe - B^We = F 
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(10.4) 



follows since all the operators arc bounded on H'^{^1), sec (3.15), so its just an ordinary differential 
equation in H'^{Q), for any r > 0. Additional regularity in time follows by applying time derivatives. 
This was proven in [LI]. Lowering the indices in (10.5): 

(10.6) GWe + A^We - BoWe - B^We = GF 

Let >C^, I E T, stand for a product of modified Lie derivatives, see section 4, of |/| vector fields in T 
and let Wei = C^We- If we repeatedly apply Lie derivatives jCt and the projection, see section 4, 

(10.7) c/^ {Gj^ W,i, + A]^ Wei, - B^i^ Wei, B,j^ Wei, - Gj^ Fi,)=0 

where the sum is over all combination of Ii + 12 = I and c/^^^ = 1. If we raise the indices again we get 

(10.8) Wei + A'Wei = -c/^'^ [Gi.Wel, + A'j^Wi,) + c/^'^ [Boi.Wel, + Bu.Wel, + Gi.Fi,) 

where c/^"^^ = 1, if II2I < |/|, and c/^^^ = if I/2I = |/|. 
Let us define energies 

(10.9) E,=E{Wel) = {WeI,Wel) + {WeI,{A' +I)Wel), EJ = 

\I\<s,IeT 

Note that in the sum we also included all time derivatives Cot- The reason for this is that when 
calculating commutators second order time derivatives show up in the first term on the right in (10.7). 
As for (3.38) we get by differentiating (10.9) 

(10.10) Ej = 2{Wel, Wei + {A' + I)Wel) + {Wei, GWel) + {Wei, {A' + G)Wel) 

Now, G is a bounded operator by (3.17). The last term can be bounded by {Wi, {A'^ + I)Wei) using 
(4.43) which also holds for A^ by (4.37), and (3.13). Therefore, the last two terms are bounded by E^ , 
where r = \I\. Using (10.8) to estimate the first term we see that the norm of the last term on 
the right of (10.8) is bounded by a constant times Ej^ plus where \\F\\J = X]|j|<^ j^q- \\jO.^F\\, 

and ||F|| = {F,F)^^'^. The same is true with the first part of the first term in on the right in (10.8) 
since I/2I < |/| there and since we have included all time derivatives in the definition of EJ. It only 
remains to deal with the second part of the first term on the right of (10.8). This term comes from the 
commutators of and A^ and we add an additional term to the energy in order to pick it up. Let 

(10.11) Dj = 2cj''^Wel,AlWel,) 

where the sum is over all Ii + I2 = I, I/2I < 1-^1 and c/^^^ = 1. This term is lower order, it is again 
bounded by using (3.13) by the energies CEJ EJ_^. Furthermore 

(10.12) bj = 2c}''' {Wel,Aj^ Wei,) + {Wel,A% Wei,) + {Wel,A'j^ Wei,) 

where, by (3.13) the second to last term is bounded by GEJeJ_^ and the last term is bounded by 
CEJ^EJ, since we have included all time derivatives in the definition (10.9). Hence, we have proven 
that 

(10.13) \Ej + bj\<CE^{E^ + \\F\\J), \Dj\<CEJeJ_^, r > 0, eZ^=Q. 

Using induction and a Gronwall type of argument, see [LI], it now follows that: 
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Lemma 10.2. There is a constant C depending only on t and on x{t,y) but not on s such that 
(10.14) eJ <C f \\F\\JdT 



In fact, integrating the first inequality in (10.13) from to t, using that Ei{0) = -0/(0) = 0, summing 
over |/| < r, and using the second inequahty gives (-B^)^ < CEJ EJ_-^+C E'^ (E'^ + \\F\\J) dr. Hence 

(10.15) Er<CEr-i + c[ (Er + \\F\\J)dT, where Er{t) = supErir). 

Jo 0<T<t 

Introducing Mr = Jq E^ dr, gives Mr — CMr < CEr-i +C c/r. Multiplying by the integrating 

factor e~*-^*, we see that Mr is bounded by some constant depending on t times CEr-i + C \\F\\J dr. 
Hence for some other constant Er < CEr-i +C Jg ||F||^dr and (10.14) follows by induction. 

From the uniform energy bounds in Lemma 10.2 it follows that \\Ws\\ < C, where C is independent 
of £ so we can choose a weakly convergence subsequence We^ that converges weakly in the inner product 
to W which is also in that space. Let U he a smooth divergence free vector field which < t < T in 
the support. Then 

(10.16) / / gab{L\W^)U''KdydT= [ [ gabW^{Ll*U') ndy dr 
Jo Jn Jo Jn 

where Lf* is the space time adjoint. The only term that depends on e in Lf* is A^, since is self 
adjoint. Since the projection is continuous it also follows that A^U — AU, as £ — ^ 0, strongly in if 
U e H^. Then right hand side of (10.16) therefore convergences so we get 

rp rj-i 

(10.17) / / gabW''iLlU'')KdydT= [ [ gabF'^U'' ndy dr 

Jo Jn Jo Jn 

where now W is the weak limit. Hence I^ is a weak solution of the equation. Furthermore is 
divergence free so it follows that W is weakly divergence free, i.e. 

(10.18) / / W''idaq)KdydT = 

Jo Jn 

for all functions smooth functions q that vanishes on the boundary. We now want to conclude that W 
has additional regularity so we can integrate by parts back and conclude that is a regular solution. 

Note that since the curl of a gradient vanishes 

(10.19) cmlA'^We = 0, when d{y) > e, 

It follows that the formulas in Lemma 9.1 hold true for d{y) > e and hence 
Lemma 10.3. When d{y) > £ we have 

(10.20) |Acurlu;e|^_i < C{\W,\^ + liy^l") + I curlF|^_i 

(10.21) \Dtcmlwe\^_^ < C{\We\'^! + \We\^) 



and by Lemma 5.2, see the last statement: 
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Lemma 10.4. 



(10.22) 
(10.23) 



<C{\ cutIw\^_^ + I divVF|^_i + \W\J) 
\W\l^ < C(|curlw;|^_i + |divW^|^_i + \W\'^) 



Let C, 



iU,e 




and for r > 1 let 



(10.24) 



U,e _ 



CUrlu;e||2^r-i(j7^) + II curiae ||wr-i(f7^) 



where ||/3|| 




and = {y G ri; (i(?/) > e}. Since divV^ = divV^ = and since d{y) > e and in the domain of 
integration in (10.24) it therefore follows from Lemma 10.3 and Lemma 10.4 that 



(10.25) 



where C depends on t and x(t, y) but is independent of e. This together with Lemma 10.2 and Lemma 
10.4 now gives us uniform bounds: 

Lemma 10.5. 



We can hence pass to the limit and conclude that the limit W also satisfies the same estimate and 
therefore if we integrate by parts in (10.17) and (10.18) conclude that W in fact is a smooth solution: 

Proposition 10.6. Suppose x{t,y) is smooth and (2.7) and (2.8) hold for <t<T. Suppose also 
that F is smooth for < t < T, divF = and F vanishes to all orders as t 0. Then the modified 
linearized equation (10.1) with vanishing initial data, Wq = Wi = 0, have a smooth solution W for 
<t <T, satisfying divW^ = 0. Furthermore, the solution satisfies the estimate: 



(10.26) 




(10.27) 




r > 0. 



11. Estimates for the inverse of the modified 



LINEARIZED OPERATOR IN THE DIVERGENCE FREE CLASS 



We will now give improved estimates for the modified linearized equation 



(11.1) 



L^W = W + AW - BoW - B{W = F, 



within the divergence free class. We have 
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Theorem 11.1. Suppose that x and p are smooth, p\qq = and that the coordinate and physical 
conditions (2.8) and (2.7) hold, for < t < T . Let Li be defined by (2.4-9) and suppose that W a,nd F 
are smooth and divergence free satisfying (11-1) for <t <T. Then ifW\^_^^ = W\^_^^ = we have 

(11.2) \\W\\r + \\W\\r<K3e^'^^+''o'^^y2^r-i-s W^Ws dr, 

s=0 

s=0 •'^ s=0 



for < t <T. If in addition = then for r > 1 and < t < T we have 

(11.4) 



\\W\\r-l + \\W\\r-l + ||W^||r-l + Co||W^||r < i^4e^^(^+"° I^Iir-l-^ / (II^IU + II^IU + l|curlF||,) dr 

s=0 

Here cq > is the constant in ( 2. 7), where 
(11.5) rig = sup ns{t), 

0<t<T 

(11-6) ns{t) = \\x{t, OIU+s.oo + ||a;(t, Olls+s.oo + P(i, Olb+s.oo 

and K3 is a constant, which depends on n_i + ci, where ci is the constant in (2.8). 

For r = (11-2) is the basic energy estimate in section 3. For r > 1 it follows from first applying 
Lie derivatives with respect to space tangential vector fields to the equation and estimating the energy 
for these as well as using the evolution equation for the curl and the fact that we can estimate any 
derivative by the curl, the divergence and tangential derivatives. The difference between (11.2) and 
(10.27) is, apart from that we keep track of how the constants depend on the solution we linearize 
around, that we only have space derivatives in the norms in (11.2). The commutators in the energy 
estimate are now estimated using the curl as well as the energies of tangential derivatives. (11.3) follows 
from (11.2) using (11.1) to estimate W. (11.2) and (11-3) follows from Proposition 11.4 below. 

The importance of (11.4) is that one gets control of an additional space derivative Co||M^||r by only 
controlling an additional time derivative and the curl of the right hand side. (11.4) without the term 
co||W^||r in the left and ||curlF||s in the right, in principle follows from (11.2) applied to the equation 
one gets for W by commuting Df through Li in (11.1). The commutator term AW can in principle be 
controlled by the energy of the same order but in order not to get constants depending on A we will 
bound it using an additional space derivative. co||W^||r can be controlled as follows. Using the estimate 
without this term in (11.1) gives control of By Lemma 5.4 this gives us control of co||T^||r 

if we also control || curlt/;||r-i. We then use that there is an improved evolution equation for the curl 
which only requires control of and not ||VF||r, by Lemma 9.2. (11.4) follows from Proposition 

11.9. 

Let us rewrite (11.1) slightly 
(11.7) W + AW = H, where H = BqW + B^W + F 

and by (2.51)- (2.52) the operators Byi and Bq are on divergence free vector fields 



(11.8) = -P(5"''(A56c - iObc - 2agbc)W'), BoW^ = P {g'^' a {Dt g^c - iObc - crgbc)W^) 
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It follows from (4.47) and (4.49) that 

(11.9) Wi + AWi = Hi + Ki, Ki = Gj^'^Ai^ Wi, , Hj = Gj^'^ {Boi, Wi, + Bu, Wj,) + Fj 
where Wj = £3^, Fj = C^^F, and Ai and are given by (4.41) and (4.43). Let 

(11.10) Ei = E{Wi) = {Wi, Wi) + {Wi, {A + I)Wi) 
Then 

(11.11) Ei = 2{Wi, Wi + {A + I)Wi) + {Wi, GWi) + {Wi, (A + G)Wi) 

= 2{Wi, Ki + Hi) + {Wi, Wi) + {Wi,GWi) + {Wi, {A + G)Wi) 

The last there terms can be estimated by Ei, by (3.42), so we get 

(11.12) \Ei\ < 2^/e'i\\Ki + Hi\\ + ^3(1 + Co^)Ei 
However, in order to estimate the first term we must estimate ||-f^/|| + 

Lemma 11.2. Let Ci, Ki, for i = 1,2, 3, nis and vris he as in Definitions 5.2 and 7.1. We have 
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(11.13) \\Gl^^^W\\<Kims\\W\\, s = |/| - - |I| 

(11.14) \\Bii,W\\ < K2ms\\W\\, s = \h\, i = 0, 1 

(11.15) \\Ai^W\\<K^(m,+^\\W\\+m,\\W\\i), s = \h\ 

and if r = \I\ then 

r 

(11.16) \\Ki\\ < K3^mr+i-s\\W\\, 

s=0 

r 

(11.17) \\Hi\\<K2Y,^r-si\\W\U + \\W\\,) + \\F\\r 

s=0 

Proof The proof of (11.14) and (11.15) uses (4.41) and (4.43) for Ai and Bi, the bounds (3.9) and 
(3.3) and (7.10)-(7.11) to estimate the pressure in terms of the coordinate. The proof (11.13) also uses 
the interpolation inequalities in Lemma 6.2. (11.16) and (11.17) is just a combination of (11.14) and 
(11.15) with (11.13) and the interpolation inequalities in Lemma 6.2. Note the remark after Definition 
5.2 that \\Wi\\ < Ki{\\W\\s + |blUI|VF||) if s = A remark about the estimate (11.15) for Ai^ 
is required. By (4.43) Ai = A^j^ which can be estimated by (3.9) if we control ||Viv<S'S'^p||£,oo(gn). 
In (11.17) we claim that this will involve at most s + 2 space derivatives of the metric. In fact, 
S^p = S^p + C^^-^'^{S^'a) ■ ■ ■ {S^''-^a)S^''p) and S^'^p = on the boundary so it follows that the 
normal derivative must fall on S^'^p so the factors S^^a never gets differentiated by Vjv- D 



Definition 11.1.. Let 



(11.18) Ef = { Ei) : = ||curlt(;||^_i + || curM||^_i, {W)A,r = J2 iWi,AWi) 

|/|<r-,Se5 \I\<r,IeS 

where G^ should be interpreted as 0. 



Summing up the results in Lemma 11.2, Lemma 9.2 and Lemma 5.3 we have: 
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Lemma 11.3. We have 

r 

(11.19) \Ef\ < Ks{l + c^^)Ef + {K2mr-s\\W\\s + K^mr+i-s\\W\l) + 



8 = 

and 



(11.20) \C'^\ + \\cni\^\r-i < K2Ymr-s{\\M\s + W\\s) + K^Ymr-s\\F[ 

s=0 s=0 

and 

r 

(11.21) lll^ll, + \\W\\r + {W)A,r <K^Y1 ^r-s{C^ + ) 

s=0 



(11.22) C:!^ + Ef <KiY,^r-s{\\W\U + \\Wl) + {W)A,s 

s=0 

Proof. The first inequality follows from (11.012) and Lemma 11.2, the second from Lemma 9.2 and the 
third from Lemma 5.3. The last inequality is just due to that contains ||KVF||f and differentiating 
K produces lower order terms. □ 

Proposition 11.4. Let cq > and ci < oo be such that (2.7) and (2.8) hold and x is smooth for 
<t <T. Let rhg = supo<t<r ms(^); where rhs is as in Definition 6.2, and set 

(11.23) Er = \\W\\r + \\W\\r + {W)A,r 

Then for r > 0, there is K^, as in Definition 6.1, such that, for <t <T, 

(11.24) Er{t) < i^3e^^(^+'=o')^ j2^r+i-s{Es{0) + f midr), 

s=0 

and for r > 1 

r „t r — 1 

(11.25) <iC3e^^(^+=o')^(5^m,+i_,(^,(0)+ / c?r) + ^ . 

s=0 •'^ s=0 

Proof. We will prove the estimate for Er = Ef + C^, in place of E^, and in view of Lemma 11.3 
and interpolation, rhrms < Kiuir+s, the estimate for Ej. follows from this. By Lemma 11.3 and 
interpolation, rhrms < Kiihr+s, we also have 

dE "^'^ 

(11.26) < Ksil + Co^)Er + KsJ2^r+i-sEs + KiJ2^r-s\\F\\s 

s=0 s=0 
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where we also used that rhi < C3. Let a = Ks{l + Cq^). Multiplying by the integrating factor we get 



r-l 



(11.27) {Ere-"')' < e-'''K3{J2^r+i-s{Es + + mo||F||,) , 

s=0 

where this is to be interpreted as that the sum is absent if r = 0. Integration of this gives that 

(11.28) Er{t)<K3e°''^[Er{Q)+ j ( ^ m^+i_,(4 + ||i^lU) + no||F||r-) c?t) , t<T, 

s=0 

where the sum is to be interpreted as absent if r = 0. The proof of the estimate with in place of Ej. 
is by induction. Since the sum is absent if r = it follows for r = 0. In general we use the interpolation: 

rrir+i-smsJ^x-t < Crhinir+i-t < K^rhr+i-t- 

To prove the estimate for we note that by Lemma 5.3 it is bounded by the curl and the 

tangential components: 

r 

(11.29) <iCi^m,_«(||curl«}||,_i+ ^ \\Wi\\) 

s=0 |/|=s,7e5 

where the curl was estimated in Lemma 11.3 and the tangential components can be estimated using 
the equation Wj = AWj + Kj + Hj and Lemma 11.2: 

7' 

(11.30) W^lW ^ {K2mr-s\\W\\, + K3mr+l-s\\W\\s) + mo||^||r+l + lli^llr 
|/|<r,/e5 s=0 

Hence by (11.29), (11.20) and (11.30) 

r r + 1 r 

(11.31) \\W\\r < K2Y,^r-s\\W\\s + K3Y,f^r+l-s\\W\\s + KiJ2'>^r-s\\F\\s 

s=0 s=0 s=0 

(11.25) follows from this with r replaced by r — 1. □ 

We now want to get estimates for an additional time derivatives by differentiating the equation. 
This gives an estimate for the normal operator through the equation and this together with estimates 
for the curl gives the estimate for an additional space derivative that we are seeking. Recall that 

(11.32) W + AW = H, where H = B{W,W) + F. 
where B given by (2.49) or (2.63) is 

(11.33) B^{W, W) = -{gac - UJac - (7gac)W'' + a{gac - UJac)W' - daQo 

where gab = D^gab- In order to differentiate with respect to time let us now write this in the form 
GW + AW = H: 

(11.34) gacW'' + A^W = BSW, W) + g^cF"^ 
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Differentiating tfie equation gives 

(11.35) gacW" + A,W + A,W = A B,{W, W) - gacW^ + QacF" + QacF" 

We have 

(11.36) A W) = - {pac - t^ac - (^Qac) - {gac - OJac + ^(^ab - ^(jac - '(JQac) W'' 

In conclusion we get 

(11.37) W + AW + AW = Hi wliere Hi = BgW + BsW + B7W + GF + F 
where 

(11-38) BgW'' = -P{g''-{2gac - Uac - ^gac)W'') 

(11.39) B^W'' = P{g''" {2&{gac - a;„e) - 9ac + <^ac + ^9ac)W'') 

(11.40) B^W'' = P{g''' {a{gae - CJac) + K9ac - <^ac + O^ao)) W^) 

Applying vector fields to (11.37) gives 

(11.41) Wi + AWi + AWi = Hu + Ku, where Ku = -G^'^^ {Aj^Wi^ + Aj^Wi^) 
and 

(11.42) Hu = Gj^'^ {Bei,Wi,+ Bri,Wi, + Bsi.Wi,) + Fj + Gj^'^Gi^Fj, 
Let 

(11.43) Eu = E{Wi) = {Wi, Wi) + {Wi, {A + I)Wi) 
Then 

(11.44) Eu = 2{Wi, Wi + {A + I)Wi) + {Wi, GWi) + {Wj, (A + G)Wi) 

= -2{Wi, AWi) + 2{Wi, Ku + Hu) + {Wi, Wi) + {Wi, GWi) + {Wi, (i + G)Wi) 

The last three terms are estimated by Eu and the second term is estimated as before by lower energies: 

(11.45) \Eu\ < 2VEu\\AWi\\ + 2^/e^\\Ku + Hu\\ + Ksil + Co^)Ei 

However, the estimate of the first term —2{Wi,AWi) requires a couple of new observations. This term 
could be absorbed by adding 2{Wi,AWi) to the energy, which instead would produce 2{Wi, AWi) 
and 2{Wi,AWi). However, we want to have an estimate that only requires two time derivatives of 
the coordinate and this would require an estimate for A, which requires three time derivatives of the 
coordinates. Instead we will use that by Lemma 5.4 and Lemma 5.5 

(11.46) \\AWi\\ < K^iWdWiW + \\Wi\\) < K3{1 + Co^){\\AWi\\ + HcurlW/H + \\Wi\\) 

Then there appears to be a loss of regularity, but remember that we now have an estimate also for 

in the energy and by the equation (11.2) we can estimate H^VF/H < ||W7|| + \\Ki\\ + H-ff/H, where 
the last two terms were estimated in Lemma 11.2. curl is by (5.22) up to terms of lower order equal 
to C^gcmlw. At this point we have to use that we have an improved evolution equation for the curl. 
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Lemma 11.5. We have 

(11.47) \\Bu,W\\ < Ksm,\\W\\, s = \h\, z = 7,8,9, 

(11.48) \\Ai,W\\<Ks{ms+i\\W\\+m,\\W\\,), s = \h\ 

and if r = \I\ then 
(11.49) 

r 

\\Ku\\ < Ks^rhr+i-siWWlU + \\W\U) 

r r 
(11.50) \\Hu\\ < K3j2^r-s{\\W\U + \\W\U + \\W\U) + \\F\\r + K2j2^r-s\\F\\ 



Definition 11.2. Let 

(11.51) El, = { Yl ^^.0'^'' C5 = ||curH|^ + ||curl«;||^ 

|/|<r,Se5 

where w is as in Lemma 9.1. 

Summing up the results in Lemma 11.5, Lemma 9.2 and Lemma 5.5 we have: 
Lemma 11.6. We have 

(11.52) E^^,<Ks{l + c^^)E^^, 

r r+l r 

+ K-,Y,^r-s{\\W\U + \\W\U) + KsY,^r+l-s\\W\U + \\F\\r + K2Y,^r-s\\F\\ 



s=0 s=0 s=0 



r+l 



and 

(11.53) < C^^ + K2J2^r+l-s\\W\U + K,\\cm\F\\r 

s=0 

and 

r r r 

(11.54) \\W\\r < Kl J2 ^r-sEl^ + rhr-,{\\Wl + \\W\U) + K^Ymr-s\\F\U 

s=0 s=0 s=l 

and 

r 

(11.55) co||W^||.+i < i^3(C5 + Ef^, + K-,Y,^r+i-s{\\W\l + \\W\l)) 



s=0 



Proof. The energy estimate (11.52) follows from the energy estimate (11.47c) and the estimates in 
Lemma 11.5. The estimate for the curl (11.53) follows from Lemma 9.2. The estimate for W (11.54) 
uses Lemma 5.3: 

r 

(11.56) \\W\\r < K,J2mr-s{\\cuvlw\U.i+El^)+mrEl, 

s=l 
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where the estimate for the curl follows from Lemma 11.3, and we use interpolation, msihr < K^rhs+r- 
Let us now prove the additional space regularity (11.56). We have from the equation, AWi = —Wj + 
Hj + X/, and Lemma 11.2 

r 

(11.57) \\W\\f^A<El^ + K2Y,rnr^i-s{\\W\\s + \\W\\:) + \\F\\r, ||W^||f,A = E W^^s^W 
and (11.56) follows since by Lemma 5.5: 

r 

(11.58) C^\\W\\r+l<K^{C^^ + \\W\\f^A + Y,mr+l-s\\W\\s) □ 

s=0 

Proposition 11.7. Let co > and ci < oo be such that (2.7) and (2.8) hold and x is smooth for 

< t < T. Let rhg = suPq<j<j. rhs^t), where rhs is as in Definition 6.2, and set 

(11.59) Er^i = \\W\\r + llM^llr + {W) A,r + \\W\\r + {W) A,r + || curltD||r + || curlit;||r + Co||VF||r+l. 
where w is as in Lemma 9.1. Then for r > there is K4, as in Definition 6.1, such that, for <t <T, 

(11.60) ErAt) < K4e^^('+^°')'^ ^m,+i_,(E,,i(0) + / (||F|U + \\F\U + ||curlF|U) dr) 
and for r > 1 

(11.61) \\w\\r-i < K^e'^'^'+^o"^^ J2^r+i-s{Es,i{0) + I (||F|U + ||F||, + ||curlF||,)dT) 

s=0 •'^ 

+^^e^.(i+i/co)T Y,m,-,-sm\s + iii^iu). 

s=0 

Proof. The proof would be the same as the proof of Proposition 11.4 apart from that we must worry 
more about the possibility of the constant cq being small. As in the proof of Proposition 11.4 the 
estimate for £^^,1 would follow from the same estimate for Er^i = Ef-^^+Cl^^^ + Erj where Ej. = E^ + 
is as in the proof of Proposition 11.4. The critical term with cq is by Lemma 11.6 and Lemma 11.3 
bounded by the other terms plus lower order terms. Note that by Lemma 11.3 and interpolation 
X]s=o ^r+i-s-E^r bounds the lower order terms with ||Vt^||s and ||VF||s, for s < r in Lemma 11.6. By 
Lemma 11.6 and the proof of Proposition 11.4 we have 

dE '^"^ 

(11.62) — ^ < ^4(1 + Co^)Er,i + Kail + Cq ^) ^ m^+i_,^,,i + i^i ^ m^_,||F||, + C||F||^ 

where we also used that rhi < 04. Let a = ^4(1 + Cq ^). Multiplying by the integrating factor we get 

r— 1 r 

(11.63) (^r,ie""*)' < e-'''K4({l + Co')^rhr+i-sEs,i + Y,mr-s\\F\\s + 

s=0 s=0 

where this is to be interpreted as that the sum is absent if r = 0. Integration of this gives that 

nt r-1 r 

(11.64) Er,i{t) < K^e"^ (^Er,i{0) + / {{1 + Co^)Y,rhr+i-sEs,i +Y,^r-s\\F\\s + \\F\\r) dr) , 

•^0 s=Q s=0 

for t < T, where the sum is to be interpreted as absent if r = 0. The proof of the estimate (11.60) 
with Er^i in place of Er^i follows by induction from (11.64). Since the sum in (11.64) is absent if r = 
it follows that it is true for r = 0. In general we use interpolation, m^+i-s^i^s+i-t ^ Crhirhr+i-t < 
K4mr+i-t. (11.61) follows as in the proof of (11.25). □ 
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12. Existence and estimates for the inverse of the 
modified linearized operator for general vector fields 

In this section we will show that the modified linearized operator can be solved for general vector 
fields outside the divergence free class, i.e. we solve 

(12.1) L,W = F, Wl^^ = Wl^^ = 

when F is not necessarily divergence free. Below we give estimates for the solution of (12.1) that are 
good enough that the linearized operator can be considered as a lower order modification of (12.1); In 
the next section we will use these to prove existence and estimates also for the inverse of the linearized 
operator by iteration. One gets a new iterate by substituting the previous iterate into the right hand 
side of (12.3) and solving for the new iterate in the left hand side. We want estimates that are good 
enough that wc get the same regularity for the new iterate so we need estimates for (12.1) that do not 
loose regularity going from F to W. We have: 

Theorem 12.1. Let < T < cq < 1 and < ci < oo 6e such that (2. 7) -(2. 8) hold and x is smooth for 
<t <T. LetUg = supo<t<r ns(i), where Ug is as in Definition 6.2. Then the equation (12.1), with F 
smooth, has a smooth solution W, for < t < T. Furthermore, there is K4 as in Definition 6.1, such 
that, forO<t<T, 

r .t 

(12.2) \\W\\r-l + \\W\\r<Kiy2^r-s I ^^'^^ r>l 

and 

(12.3) \\W\\r-i < K4^n^_g \\F\\sdT + KA^n^_^_g\\F\\g, r>l 

s=l •'^ s=0 

As in section 3 we can decompose W = Wq + Wi where Wq is divergence free and W\ is the gradient 
of a function vanishing at the boundary. By (3.26) Wq satisfies 

(12.4) LiWo = ~AWi + B^^Wi + Bq^Wi + PF, Wo|,^o = W^o|,^o = ^ 

where all the terms in the right hand side are divergence free and -Bqi and Bn are bounded operators 
given by (3.25). By (3.27)-(3.28) Wi satisfies 

(12.5) M^i" = 5"'5bgi, Agi = v?, gi|aj, = 0, 
where 

(12.6) + a^ = divF, ^\^^^ = A<^|,=o = « 

The solution of (12.6) is a smooth function if F is smooth so it follows that that Wi is smooth and 
hence (12.4) has a smooth solution Wo by Theorem 10.1. Therefore, we have proven that the modified 
linearized operator (12.1) has a smooth solution W ii F and x are smooth and the coordinate and 
physical conditions are satisfied for < t < T. However, in the right hand side of (12.4) the term AWi 
looses space regularity since A is order one. If we just use Proposition 11.4 and Proposition 12.3 below 
we are going to get an estimate that looses space regularity going from F to in (12.1). However, 
because the curl of AWi vanishes we can use the improved estimate in Proposition 11.7 that gains an 
extra space derivative to handle the term —AWi. Let us first prove the estimate for (12.5)-(12.6): 
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Lemma 12.2. Suppose that 

(12.7) D'^ip + a<f = b^ip-2&bt(p + &'^cp = f, 

Let T < 1 and set fh^ = supg^^^rp ms{t) , where rhg is as in Definition 6.2. Then, there is K^, as in 
Definition 6.1, such that, for < t <T and r > 1 we have 

(12.8) Mr-i + M\r-i<Ksy2m,_^_, WfWsdT, 

(12.9) < £ (ll/IU + /Il/IU dr) 

Proof. (12.4) is just an ode for each space coordinate, however one just has to make sure to integrate 
it up in such a way that we do not get more than two time derivatives on the metric. 

(12.10) A {{bt<pf + CT^V^) = 2a{Dt<pf + 2a{a - a^)<p^ + 2( A</')/, AV - 2^A<^ + f> V = / 

Integrating this in time and space gives the lowest order energy estimate in (12.8) and the lowest order 
estimate in (12.9) follows from this since once we have estimates for the cp and Dttp we get an estimate 
for D^(p from the equation. 

In order to get (12.8) and (12.9) for higher derivatives we commute through , defined in section 
4 by R^ f = R^ (nf), where / = [ii, ...,ir) is a multiindcx and R^ = Ri^ • • • Ri^ is a product of the 
vector fields in TZ defined in section 4. Then [A, i?^] = and with ifj = W Lp and d/ = W a, we obtain 

(12.11) AV/ - 2abt<pi + dV/ = //, // = 2c^'^'&i,btipi, - d^^^'^'&i.&i.ifi, + R^ f 

were the sums are over all combinations of Ji + /2 = / respectively Iq + Ii + 12 = I and c^^^^ = i an^j 
(II0I1I2 — I unless I2 = I in which case they are 0. We can now use (12.10) applied to // in place of / 
and cpi in place of (p. Here the terms in // are lower order. □ 

Once we get the corresponding bounds for ip in terms of divF, the bounds for Wi follows from 
Proposition 6.1: 

Proposition 12.3. Suppose that = g^-^d^q, where q\q^ = and Aq = ip where ip satisfies 

(12.12) AV + f^V? = divF, 

Let T < 1 and set = supQ<«2^ ms(t), where rhs is as in Definition 6.2. Then, there is K^, as in 
Definition 6.1, such that, for < t <T, 

(12.13) ||I^i(t)||, + ||P^i(t)||, <iC3X^m.-.(l|^i(0)|U + ||I^i(0)||,+ / ||F||,dr), r>l 

s=l -^0 

and 

(12.14) \\Wi{t)\\r < KsTirkr-sillWimU + \\WimU + f \\F\UdT+\\F{t)\U), r>l 

s=l -^0 
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We will now further decompose the solution of (12.4) into two parts Wq = Wqq + Wqi, where 

(12.15) LiWoi = -AWi = Foi, = W^o4=o = 
and 

(12.16) LiT^oo = PF + 5nWi + SqiWi = Fqo, W^oo|,=o = ^ooL=o = ^ 

For (12.15) we use the estimate in Proposition 11.7 and for (12.16) wc use Proposition 11.4. This 
together with Proposition 12.3 gives Corollary 12.4 below. Our solution to (12.1) is now obtained as a 
sum oiW = W\ + Wqi + Woo so it will satisfy the worst of the estimates in Corollary 12.4 and this 
proves Theorem 12.1. 

Corollary 12.4. Let < T < cq < 1 and c\ < oo he such that (2.7) and (2.8) hold and x is smooth 

for < t < T . Let = supQ<j<2^ ^^(t), where rig is as in Definition 6.2. Let Wi be the solution of 
(12.5)-(12.6), let Wqi be the solution of (12.15) and let Wqq be the solutions of (12.16). Then there is 
K4, as in Definition 6.1, such that, for < t <T and r > 1, 

\\W4r + \\W4r < K,Y,Rr-l-sf\\F\\sdT, 

<K^y^n,_^_lf\\F\\,dT+\\F\l) 



(12.17) 



\W- 



1 r 



and 



(12.18) IIT^oi 11.-1 + II Woi||r-i + 11^01 II. <i^4Vn,_, / \\F\ldT 



and 



(12.19) 



1 

|W"oo||r + llW^Oollr < ^4 / WHsdr, 

r ,.t r-1 

||Woo||r-l < i^4^n^_l_« / \\F\\,dT + K^Y,^r-l-s\W\\ 
s=0 •'^ s=0 



Proof. (12.17) follows from Proposition 12.3. By Proposition 11.7 we have for r > 1 
(12.20) 

||W^Ol||r-l + ||W^Ol||r-l +C0IIW01II.+ ||Woi||._l < -fC4 Vzir-l-s / (l|-^Ol|U + l|-^Ol|U + II CUrlFoi || «) cZt 

s=0 



We remark that it follows that also Woi|j=o = since = 0. Here the curl of Fqi = AWi 

vanishes and DtAWi = AWi + AWi - GAWi so 

r 

(12-21) l|i^Ol||r-l + ||i^Ol||r-l < J]n._i_«(||Wi|U + llWilU) 

s=0 
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Using (12.17), (12.20) and (12.21) we obtain (12.18). Note that the constant cq in (12.20) can be 
replaced by 1 since we have two consecutive integrals and we assumed that < t < T < cq. Finally 
from Proposition 11.4 we get 

sdr, 

s cLt + i^3||Foo||r-l 

Now, the operators Bqi and Bn in (12.16) are bounded operator given by (3.25), of the same form that 
we have already studied in section 9 and PF, the projection is bounded by the estimates in Proposition 
6.38, so it follows that 

T r 

(12.23) llFoollr < K^^Tir-siWWiW, + llVFill,) + i^i^n^_,_2||F||, 

s=0 s=0 

Combining these inequalities, using interpolation as usual, gives also (12.19). □ 

13. Existence and estimates for the inverse of the linearized operator. 
In this section we finally prove existence and estimates for the the inverse of the linearized operator 

(13.1) LoW = F Wl^^ = Wl^^ = 0, 
where Lq = ^'{x) is given by (2.14). (13.1) can be written 

(13.2) LiW = B3W + F W^^^ = W\^^^ = 

where the modified linearized operator Li is given by (2.49) and is given by (2.57). In the previous 
section we proved existence and estimates for the modified linearized operator Li: 

(13.3) L,W = F wl^^ = Wl^^=0, 

The existence and estimates for (13.3) can now be used to prove existence and estimate for (13.2), and 
hence for (13.1), by iteration. We simple define a sequence by Wq = and for A; > 1: 

(13.4) L,Wk = BsWk-i + F Wkl^^ = Wkl^^ = 

We will use the estimates for (13.3) to show that Wk converges to a solution of (13.2) and that the 
solution of (13.2) satisfies the same as estimates as the solution of (13.3). 

Theorem 13.1. Let < T < cq < 1 and < ci < oo 6e such that (2.7)- (2. 8) hold and x is smooth for 
<t <T. Let rig = supo<t<r ns(i), where rig is as in Definition 6.2. Then the equation (13.1), with F 
smooth, has a smooth solution W, for < t <T. Furthermore, there is K4 as in Definition 6.1, such 
that, forO<t<T, 

(13.5) \\W\\r-i + \\W\\r<K4i2l^r-s [ \\F\\sdT, r>l 
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OOlk "f II ►'•'OOllr ^ -'»-3 



(12.22) 



00 



nt 

r-1 < Vnr-l-s / ll-^C 
.-n •''0 



and 

r-l 



(13.6) 



\W\\r-i<Ki^n^_l \\F\\,dT + K4j2nr-i-sn, r>l 



Proof. The existence and the estimates in the theorem for (13.3) were given in Theorem 12.1. The 
estimate for (13.1) follows from the estimate for (13.3) by writing (13.1) in the form (13.2). If W 
satisfies 

(13.7) LiW = BsW + F, wl^^^ = Wl^^ = 

where B3 is given by (2.57), then by (13.5) for (13.3) 



(13.8) 



r-l 



+ \\W\\r<K^i2n,_, [ {\\F\\, + \\W\U)dT, r>l. 



We claim that (13.5) for (13.3) follows from this with W = W, by induction, for some other K4. In 
fact, assume that (13.5) is true foi r < k — 1, then it follows from (13.8) and interpolation that 

(13.9) \\W\\r-l + \\W\\r < K^y^^r-s [ \\F\\sdT + K4 [ \\W\\rdT 

+ Kiy^yZ^r-snLs-k / \\F\\kdz dT < K^y^nLr-s \\F\\^ dr + \\W\\rdT 
s=lk=l «=1 -^0 



for some other K/^. By a standard Gronwall argument wc can get rid of the ||l^||r; replacing by 
some other K4. Let g{t) = \\W\\r dr and f{t) = J2l=oI^-r-s Jq II-^IU ^t. Then g'{t) < K^g + K^f so 
{ge~^*^)' < K^f and integrating this up gives g < K4 f dr for some other K4 and for t <T. 
Similarly it follows from (13.5) that the solution of (13.7) satisfies 



r-l 



(13.10) \\W\\r-i < K^y^^r-s (II^IU + l|W^IU)(iT + X4^n.-i-3(l|i^lU + l|W^IU), r>l 

s=l •'^ s=0 

(13.6) for (13.1) follows from this with W = W using the estimate (13.5) that we used obtained. 

It remains to prove existence for (12.2). We put up an iteration Wq = and LiWk = F + (Li — 
Lo)Wk-i, for A; > 1. Then LiWi = F so Wi satisfies the desired estimate and is smooth. Let 
Wk = Wk- Wk-i, for A; > 1. Then Wi = Wi and LiWk = {Li - Lo)Wk-i, for A > 2. In conclusion 

(13.11) L^W,=F, L,Wk = B3Wk-u k>2, Wkl,^^ =Wk\t^o = ^ 
where B3 is a bounded operator given by (2.57). Using the estimate (13.8) for each k 

(13.12) ^^^P {0k{r, -nr-l + \\Wk{T, -nr) <K4j2llr-s (II^IU + ^ ^ ^ 
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Note that the supremum is inside the sum since we use (13.8) for each Wk and since in the left of (13.8) 
we may take the supremum of r < t. The same argument that lead to the proof of the estimate (13.5) 
for (13.1) from (13.8) now gives the uniform estimate 

(13.13) V sup (||T^fc(r,-)||r-i + ||W^fc(r,-)l|r-) <if4Vn,_, / \\F\UdT, r > 1, 

k=lO<r<t Jo 

where K4 is independent of A'^. (One replaces the sum in the right of (13.12) by the larger sum in the 
left of (13.12).) Similarly the uniform estimates corresponding to (13.5) also hold as is seen by using 
(13.10) for each k and replacing the sum in the right by the larger sum in the left and using (13.13) 

(13.14) V sup ||#fe(T,-)||r-i <i^4 Vn,_, / \\F\\,dT + K^y^n,_^_, sup ||F(t,-)|U, r>l 

t^l^<^<* 7^1 Jo t^Q 0<T<t 

It follows that Wn = J2^=i is a Cauchy sequence in C^([0, T],H'^~^ {^)), for any T, and hence there 
is a limit W G C'^{[0,T], H'^~^ {Q)), for any T. Additional regularity in time follows from differentiating 
this equation. We have already proved that D^W = AW + BqW + B^W + B3W £ C^i[0,T],H''-^(n)), 
i.e. DfW is continuously differentiable with respect to time so G (^[0,T], H^~'^ (Q,)) and so on. 
Since this argument is true for any r it follows that W is smooth. □ 

14. Estimates for the physical and coordinate conditions. 

We assume that the physical condition and the coordinate condition hold initially at time for 
some constants cq > and ci < 00 and we need to show that this implies that they will hold with cq 
replaced by cq/2 and ci replaced by 2ci, for < t < T, if T is sufficiently small. 

Let us introduce the space time norms: 

(14.1) lll^lllr = sup ||tt(i, •)||r,oo) II I'^^ll l^fc = || K|| |r + ••• + || |r 

0<t<T 

We have: 

Lemma 14.1. Let M{t) = supy^^ ^\dx/dy\^ + \dy/dx\^. Then 

(14.2) M(t)<2M(0), for t<T, if r|||i;|||iM(0) < 1/8 

Let N{t) = supy^QQ \Vnp\~^- Then assuming that T is so small that (14-2) hold we have 

(14.3) N{t)<2N{0) for t<T, if r|||p|||iM(0)7V(0) < 1/8 

Proof We have \Dtdx/dy\ < \\\x\\\i and \Dtdy/dx\ < \dy/dx\'^\Dtdx/dy\ soM'(t) < (1 + M2)|||i;|||i < 
2M2|||i;|||i, since also M{t) > 1. Hence 

(14.4) M(t) < M(0)(1 -2|||i;|||iM(0)t)"\ when 2|||x|||iM(0)t < 1. 

Now, VnP = N°-daP, where N is the unit normal, so DiVnP = V^Dtp + {DtN"-)daP = S/nDiP + 
{DtN^)gabN''VNP, since p\g^ = 0. Furthermore = A(5af,iV"iV^) = 2gab{DtN'')N^ + (Affab)A^"iV^ 
and iV" = {dy''/dx')N\ where 6ijN'N^ = 1. Hence IAVatj?] < M{\dDtp\ + \dDtx\\VNp\) Therefore if 
N{t) = supy^g^ |VivPr\ we have N' < M|||p|||iiV2 + M|||i;|||iiV/2 and if we use (14.2) and multiply 
with the integrating factor, N{t) = Ar(t)e-*^(°)lll^llli we get N' < 2e^/^M(0)|| Hence 

(14.5) N{t) < N{0)e^/^{1 - iV(0)2e^/^M(0)|||;3|||it)",' when iV(0)2e^/^M(0)|||p|||it < 1 
This proves the lemma. □ 



It now follows from Lemma 14.2: 
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Lemma 14.2. Let xq be the approximate solution satisfying (2.12) and suppose that (2.7) and (2.8) 
holds when t = 0. Then there is a Tq > 0, depending only and an upper hound for |||.xo|||4,2) ci o-nd Cq^ 
such (2.7) and (2.8) hold for < t < T with cq replaced by cq/2 and ci replaced by 2ci provided that 

(14.6) < T < To, and — a;o|||4,2 < 1, o,nd ~ ^o)\t=o~ ^ti^ ~ ^^)\t=o~ ^ 

Proof. We need to satisfy the conditions (14.2) and (14.3) in Lemma 14.1. Since |||x|||i < |||a^o|||4,2 + 1 

(14.2) hold if T < (8ci(|||a:o|||4,2 + ^))~^ ■ To satisfy (14.3) we use the estimate in Lemma 6.4, where 

is as in Definition 6.1, to obtain ||p||i,oo < -f^(||2;||3,oo + ||2;||2,oo + ||3j||i,oo) for some increasing function 
F. Hence (14.3) hold if r< co(8ci + 'f(|||xo|||4,2 + 1))"^- □ 

15. Tame estimates for the inverse of the linearized operator. 

We are now going to modify the estimate for the inverse of the linearized operator in Theorem 13.1 
so it can be used with the Nash-Moser inverse function theorem in section 18. We want tame estimates 
for the inverse of the linearized operator 

(15.1) ^'{x)5x = 5^, 0<t<T, (5x|^^Q = A Sx\^^^ = 

but the norms in Theorem 13.1 are in terms of W"' = 6x^dy°' /dx"^ and = d^^dy'^/dx^ and we like to 
see our operator as an operator on 6x. Using interpolation and Theorem 13.1 we get 

(15.2) 

+ + < K2(||T^||r+||'W^||r+||'t^||r)+i^2(||x||^+l + ||ir||^+l + ||x||^+l)(||W 

< Ki sup \\F{t, ■)\\r+l + sup (||x(r, •)||r+4,oo + ||a;(T, •)||r+4,oo + |k(T,-)||r+4,oo) SUp ||F(r,-)||i 

0<T<t 0<r<t 0<T<t 

< K4, sup \\6<^{T,-)\\r+i+K4^ sup (||x(r, •)||^+4,oo + ||a;(T, •)||r+4,oo + ||a;(r, •)||r+4,cx.) sup ||5$(t, Olli 

0<T<t 0<T<t 0<T<t 

Another issue is that we have estimates of Sx but we need L°° estimates for x. The norm is 

bounded by the L°° norm and the L°° is by Sobolcv's lemma bounded by the norm of an additional 
n/2 derivatives so one can obviously turn one into the other with an additional loss: 

(15.3) ||n(t,-)||r < C^|k(t,-)||r,oo < a||n(t,-)||r+rn, Tq = [n/2] + 1 

Furthermore, the Nash-Moser theorems that we will follow are in terms of Holder spaces, but one can 
obviously also turn Holder norms into L°° norms with a loss of an additional derivate: 

(15.4) C-^Mt, ■)\\k,oo < Mt, ■)\\a,oo < Ck\\u{t, OIU+i.oo, k<a<k + l 
where \\u(t, •)||a,oo denotes the Holder norms in section 17. Let us now introduce the norms 

(15.5) = + ... + |||£'t''M|||a,, where |||?x|||a = sup \\u{t,-)\\a,oo 

0<t<T 

It follows that if (2.7) and (2.8) hold then (15.1) has a solution that satisfies 

(15.6) lll<5^ll|a,2 < Ki{\\\Sma+ro+2 + |||'5$|||l || |x|| , a > 

We in fact want to solve for u in (2.13): 

(15.7) $(u) = $(« + xo) - $(xo) = fs 

Then ^'{u) = ^'{u + Xq) and the norm of x in (15.6) may be replaced by the norm oi u = x — Xq since 

(15-8) llkll|a,2 < Ilk - Xo|||a,2 + ||ko|||a,2 < " a;o|||a,2 + Ca 

for some constant Ca depending on xq. Hence we have proven: 
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Proposition 15.1. Suppose that x is smooth for < t < T and that the conditions in Lemma 14-2 

hold. Then if 5^ is smooth for < t < T (15.1) has a smooth solution 6x. Furthermore there are 
constants Ca, depending on the approximate solution xq, on (co,ci) in (2.7)-(2,8) and on a, such that 

(15.9) lll^a;|||„,2 < C„(|||(5$||U+^„+2 + - xo||| a+ro+6,2 
provided that 

(15.10) ||k-a;o|||4,2 < 1 

16. Regularity properties of the Euler map and tame 
estimates for the second variational derivative. 



Recall that the Euler map is given by 



(16.1) ^{x)i = D^Xi + dip, in [0,T] x ^, where di = -^da 



where p = *(x) is given by solving 

(16.2) Ap = -{diV'')dkV\ V' = Dtx\ p\ 



an 



We will now discuss the regularity properties of $ needed and the definition of derivatives of 
Let T = C°°([0,r] X n), J^M = {x G JT; \dx/dy\ + \dy/dx\ < M} and let 4 = / x • • • x / be /c copies 
of I = [— e, e], £ > 0. Suppose that x G C'^{Ik,J^M), m> k then we claim that $(x) G C"''(/fe, J"). In 
fact, by the proof of Lemma 7.3 p = '^{x) G C"^{Ik,J^), since there t e R was just any parameter and 
we can replace it by f G R'^ and replace the derivatives with respect to t by partial derivatives. 

Definition 16.1. Suppose that x ^ T = C°°([0,T] x O) and Wj G T., for j < k. Set x = x + r\W\ + 
... + rfcWfc and suppose that $(x) is a function of (ri, r^) close to (0, 0) with values in T . We 
define the fc:th (directional) derivative of $ at the point x in the directions ■iWj, i = 1, .., by 

d d 

(16.3) $('=)(a;)(u;i, ...,Wk) = ' ' ' ^^(^)|^,=...=^,=o' x = x + nwi + ... + VkWk 

We say that ^{x) is k times diffcrcntiable at x if ^(x) is a function of (ri, ...,rk) close to (0, ...,0) 
with values in J^, and if ^^^\x){wi, ...,Wj) is linear in each of the arguments wi, ...,Wj, for j < k. 

It is clear that (16.3) is independent of the order of differentiation, debut to conclude that it is multi 
linear in wi, ...,Wk one also needs to assume that it is continuous as a functional of x, wi, ...Wk, see [Ha]. 
We instead take (16.3) as the definition of the derivative and once we calculated it the linearity follows 
by inspection in our case. We will assume that $ is twice differentiable in which case it follows from 
the above definition that Taylor's formula with integral reminder of order two hold: 



(16.4) {(!>'{v) -^'{u))w = <^"{u + s{v -u)){v -u,w)ds 

Jo 

(16.5) ^{v) -^{u) -^'{u){v -u) = {l-s)^"{u + s{v-u)){v-u,v-u)ds 

Jo 

The Nash-Moser technique uses these reminder formulas together with tame estimates for the second 
variational derivative that we now will derive: 
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Proposition 16.1. Suppose that x is smooth for < t < T and that the conditions in Lemma 14-3 

hold. Then ^ is twice dijjerentiable and the second derivative satisfies the estimates 

(16.6) |||$"((5a;,ea;)|||„ < Ca (^|||fe|||a+4,i|||ex|||i,i + |||5a;|||i,i|||ex|||a+4,i) 

+ Ca|||x - Xo|||a+5,l|ll'^^ll|l,l|lk3:|||i,i^ 

provided that 

(16.7) II |x - a;o|||4,2 < 1 
Here the norms are as in (15.5). 

Let us now calculate the second derivative of $ and afterwards prove the tame estimates for it. Let 
us first recall the commutator identities: 

Lemma 16.2. We have 

(16.8) [5,di] = -{dM)dk 

(16.9) [5, didj] = -{dM)djdk - {dj5x^)didk - {didj6x'')dk 

Furthermore 

(16.10) [(5, A] = -{l\5x^)dk - 2{d'Sx^)didj 
and if e is another variation then 

(16.11) [5, [e, A]] = {{A5x^)diex'' + {didmSx'')diex"' + {l\ex^)dM + {didmex^)di5x'^)dk 

2{{d^5x"')dmex^ + {d''ex"')dmdx^ + {d"'Sx'')dmex^)dkdi 

Proof. (16.8) was proven in Lemma 2.2 and (16.9) follows from this since [6,didj] = [6,di]dj +di[d,dj]. 
(16.10) follows from contracting (16.9). (16.11) follows from using (16.9) and (16.10) applied to S as 
well as e in place of S. □ 

Let x{t, y, r) = x{t, y) +r 5x{t, y). The first variational derivative $'(x) of the Euler map 

(16.12) ^'{x)5xi = 5^{x)i = I 

or 

is given by 



lr=0 



Lemma 16.3. 

(16.13) ^'{x)6xi = D^6xi - dkP diSx'' + dip'{6x), 
Here 6p =p'{6x) = 'if'{x)6x satisfies 

(16.14) ASp = 5Ap + dkpASx''+2{didkp)d'6x'', where 

(16.15) 5Ap = 2dkV' didx^ diV''- 2dkV' diSv'' 
where 6v = D^dx and Sp\g^ = 0. 

Proof. This follows from a calculation using that S — Sx^dk commutes with di and hence with A or 
using (16.9). □ 

Let x{t, y, r, s) = x{t, y) + rSx{t, y) + s exit, y). Then the second variational derivative is given by 

(16.16) ^"{x){5x,ex)i = e5^{x)i = ^^^1 

or OS \r=s=Q 

is given by: 
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Lemma 16.4. Let Sv = D^Sx and ev = Dtex. Then 

(16.17) ^"{5x,ex)i = dkP {diex^ di6x'' +di6x^ diex'') -dkp'{ex) diSx^ - d^p' {5x) diex^ + di p" {5x , ex) 
where Sp = p'{Sx) = '^'{x)5x and Sep = p"{Sx, ex) = '^"{x){Sx, ex) satisfies 

(16.18) A{6ep) = [A, 5e]p + SeAp, [A, 5c]p = h + 2/2 - /s - 2/4, 5eAp = -2/5 + 2^ - 2/7 
where: 

h = {A5x'){d,€p) + {Aex'){d,5p) 
f2 = {d,d,Sp){d,ex') + {d,d,ep){d,6x') 
fs = djp{{di6x^){Aex') + {diex^){A6x') + 2{dk6x'){dkdiex^) + 2{dkex'){dkdi6x^)} 
U = didjp{{dkdx^){dkex') + {dkSx'){djex'') + {dkex'){djdx'')} 
h = {dkv'){div^){{dM){djex') + {diex^){dj5x')} + {div^){djv'){dk5xi){diex') 
k = {dkv'){{di5v''){d,ex') + {d,ev''){dj5x') + {d,5v'){diex'') + {djev'){di5x'')} 

/7 = {di5v^){d^ev') 

and Sep\Q^ = 0. 

Proof. A calculation using that [S,di] = —{di5x^)dk and e5x = gives (16.17). (16.18) follows from 
using Lemma 16.2 and 

(16.19) ASep = [S, A]ep + [e, A]Sp + [6, [e, A]] + SeAp □ 

The estimates for the first and second derivative of p = '^{x) are given in the following lemma: 

Lemma 16.5. Let p = ^'(x) be the solution of Ap = —{diV^)djV^ , p\g^ = 0, where V = Dtx. Let 
Sp = p'{Sx) = '^'{x)Sx be the variational derivative. We have with DfSx = Sv, D^ex = ev: 

(16.20) \\Sp\\r,^ < i^3(||5^||.,oo + \\Sx\Ul,^ + {\\X\\ 

r+2,00 

){\\Sx\\i,^ + \\Sv\\i,^)) 
and with p" {6x, ex) = {x){6x, ex) the second variational derivative, we have 

(16.21) \\p"{Sx,ex)\\r,oo < K3{\\Sv\\r+i,oo + ||5a;||r+2,oo)(||ea;||i,oo + He-vHi.oo) 

K3{\\ev\\r+l,oo + ||ea;||r+2,oo)(||'5x||i,oo + ll'J'vlll.oo) 

+ -f'^3(||^^||r+2,oo + ||a;||r+3,oo)(||ea;||i,oo + ||e^^||i,oo)(||(^a;||i,oo + ||(^^^||i,oo) 

Proof. The proof of (6.20) is similar to the estimate of a time derivative in the proof of Lemma 6.4. By 
Lemma 16.3, Lemma 6.3 and Lemma 6.4 

(16.22) 

\\A6p - SAp\\r-i,oo < ^i||'^a;||^+i,oo||p||i,oo+^i||p||r+i,cxD||5a:||i,oo + i^i||a;||r+i,oo||p||i,oo||^ 

< i^sll^xll^+i^oo + (||^^||r+l,oo + ||a;||r-+2,oo)||<^a;||i,oo 
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and 

(16.23) \\SAp\\r-l,oo<K3{\\6x\\r,o^ + \\6v\\r,oo)+K3{\\ 

v\\r,oo + \\x\\r,oo) ( || (^3^ || l,oo + ||l,oo) 

which proves (16.20). Similarly by Lemma 16.4, Lemma 6.3, Lemma 6.4 and (6.20) 

(16.24) ||[A,(5e]p||^_i,oo < -fs^i||'5a;||r+i,oo||ep||i,oo+^i||ep||r+i,oo||'^a;||i,oo+^i|k||r+i,oo||ep||i,^ 

r+l,oo ||c2;||l,oo + ^l|| r+l,oo ll^-'^ll l,oo l.oo 

||r+l,oo + Ibll \\Sx\U ,oo 

< ^1(116x11^+2,00 + ||e'y||r+i,oo)||<5a;||i,oo + -f^i(||<5a;||r+2,oo + ll^i'll r+l,oo) ll^^llljOO 

+ -f^'l(|k||r+3,oo + ||'y||r+2,oo)||ex||i,oo||<5a;||i,oo 

and 
(16.25) 

\\d€Ap\\r-l,oo < K3{\\6x\\r,oo + \\^'>^\\r,oo) ( || e2;|| 1,00 + || || l,oo)+-?^3 ( || £2; ||r,oo + || ||r,oo) (||fe||i,oo + ||5i;||l,oo) 

+ -fir3(||'y||r,oo+ Ikllr.oo) ( || 1,00 + || f^^" || l,oo) (|| £^||l,oo ~t~ ll^^llljOo) 

which proves (16.21). □ 

It now follows from Lemma 16.1, Lemma 16.6, the fact that di = {dy"' /dx'^)d/dy°' and interpolation: 
Lemma 16.6. 

(6.26) \\^"{ex,6x)i\\r,oo < K3{\\6v\\r+2,cK> + ||5a;||r+3,oo)(||ea;||i,oo + ||e^^||i,oo) 

K3{\\ev\\r+2,oo + hx\\r+3,oo)i\\Sx\\i^oo + ||^^^||l,oo) 

+ -f'^3(ll^^l|r+3,oo + ||a;||r+4,oo)(||ea;||i,oo + ||et'||i,oo)(||<^a;||i,oo + ||^^^||i,oo) 



Finally, also using (15.8) we get Proposition 16.1. 



17. The smoothing operators. 



We will work in Holder spaces since the standard proofs of the Nash-Moser theorem uses Holder 
spaces. The Holder norms for functions defined on a compact convex set B are given by, if /c < a < /c+l., 
where A; > is an integer, 

/-.^-.x II II II II \d"u(x) — d°'u(y)\ , , 

(17.1) |F||a,oo = = sup 2^ 1 _ \a-k + ^UP \u{x)\ 

and ||n||j:/o = sup^g^ l^^(^)l- Since we use the same notation for the norms, ||it||fc,oo = ll'^llc'' we will 
differ these by simply using letters a, b, c, d, e, f etc for the Holder norms and i, j, k.l, .. for the C'^ norms. 
However, since a Lipschitz continuous function is differentiable almost everywhere and the norm of the 
derivative at these points is bounded by the Lipschitz constant, we conclude that for integer values this 
is the same as the L°° norm of d"u for |a| < k, and furthermore, since all our functions are smooth it 
is the same as the supremum norm. Our tame estimates for the inverse of the linearized operator and 
the second variational derivative are only for norms with integer exponents, with B = Q. However, 
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since ||tt||fc,oo < C||'w||a,oo < C'||^t||fe+i,oo5 if A; < a < A; + 1, see (17.2), they also hold for non integer 
values with a loss of one more derivative. 

The Holder norms satisfy 

(17.2) \\u\\a,oo < C\\u\\b,oo, a<b 
and they also satisfy the interpolation inequality 

(17.3) ll^^llc,oo<C||u||^,^||n||^_;„^ 

where a<c< 6, 0<A<1 and Aa + (1 — X)b = c. 

We will use norms which consist of Holder norms in space and supremum C'' norms only in time 

(17.4) lll«ll|a,fc= sup {\\u{t,-)\\a,oc + \\Dtu{t,-)\\a,oo + ■■■ + \\D'^u{t,■)\\a,oc)■ 

0<t<T 

For the Nash-Moser technique, apart from tame estimates one also needs smoothing operator Sg 
that satisfy the properties below with respect to the Holder norms, and in fact also with respect to the 
norms above since the smoothing operators will be invariant under time translation. We have: 

Proposition 17.1. 



(17.5) ||5'0u||a,oo < C'||u||b,oo, a<b 

(17.6) \\Seu\\a,oo<Ce''-''\\u\\b,oc, a>b 

(17.7) \\{I-Se)u\\a,oo<C9''-'\\u\\b,oo, a<b 

(17.8) \\iS29 - Se)u\\a,o. < Ce''-''\\u\\b,oo a,b>Q. 



where the constants C only depend on the dimension and an upper bound for a and b. 

Moreover, these estimates hold with the norms replaced by the norms (17.4) for fixed k. 

First we note that (17.8) follows from (17.6), when a>b and (17.7), when a < b. (This alternatively 
follows from an additional property \\d/d9S0u\\a,oo < C'^"~^~^||'"||6,oo) a>0. that also hold.) 

For compactly supported functions on R" there are standard smoothing operators, sec [HI], that 
satisfy the above properties (17.5)-(17.8), with respect to the norms defined in (17.1). However we have 
functions defined on the compact set Q, that do not have compact support in Q. Therefore we need 
to extend these functions to have compact support in some larger set, without increasing the Holder 
norms more than with a multiplicative constant. There is a standard extension operator in [S] that 
turns out to have these properties, see Lemma 17.2 below. If Sq is the standard smoothing operator 
mentioned above, that satisfies (17.5)-(17.8), the we define our smoothing operator by 

(17.9) Sou = S0u\^, where u = £xt{u) 

Since Sq satisfies (17.5)-(17.8) and since HuHb^oo ^ C'||'fi||6,oo) by Lemma 17.2, it follows that Sq satisfies 
(17.5)-(17.8). 
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Lemma 17.2. There is a linear extension operator Sxt such that Sxt{f) = f in {y; \y\ < 1}, supp 
Sxt{f) C {y; \y\ < 2} and 



(17.10) \\Sxt{f)\\a,oo < C 

where the norms in the left are Holder norms in {y; \y\ < 2} and the norms in the right are Holder 
norms in {y; \y\ < 1}, and C is bounded when a is bounded. 

Proof. We will introduce polar coordinates and for fixed angular variables co extend a function defined 
for the radial variable r < 1 to r > 1. Away from the origin, the change of variables given by 
polar coordinates is a diffeomorphism and Holder continuity is preserved under composition with a 
diffeomorphism k: 

(17.11) \\f O K\\a,oo < Ca\\f\\a,oo 

Therefore, let us first remove the origin by a partition of unity. Let xo G C^(R) satisfy Xo(|2/|) = 1; 
when \y\ < 1/2 and Xo(|y|) = 0, when \y\ > 3/4, and let xi = 1 — Xo- Furthermore, we multiply with 
another cutoff function so that the extension has compact support in |y| < 2. Let X2 € C^(R) satisfy 
X2(|y|) = 1, when \y\ < 5/4 and X2(y) = 0, when \y\ > 3/2. If £xti{f) is the extension operator in the 
radial variable, defined in (17.14) below, we now define the extension £xt{f) of / to be 

(17.12) Sxtif) = X2Sxt,{xif) + xo/ 

Holder continuity in {r,u>) follows from Holder continuity of £xti{f) in the radial variable and the 

linearity and invariance under rotations of £xti{f), using the triangle inequality. In fact if fi^{r) = 
f{r,uj) then d^£xti{f^) = £xti{f2), where = d^f^ and if j + |a| = A; < a < A; + 1 then by (17.18) 
and (17.17) 

(17.13) \d:l£xt,{fS){r) - di.£xt,{f^){p)\ < xti(/«)(r) - dl£xt,{f^){p)\ + \dl£xt,{f^ - f^){p)\ 

^ ia^ag/(r»-a^ag/(p-,.;)i _ ^ \dm-iP,^')-didsfipy)\ _ . 

r',p' ^ |w'-(t'|«-'= 

It therefore remains to prove the estimates (17.17) and (17.18) for the extension in the radial variable 
only given by (17.14). 

Suppose that /(r) is a function defined for r < 1, then we define the extension / by £xti{f){r) = 
f{r), when r < 1, and 

/oo 
/(r-2A(r-l))Vi(A)cZA, r>l 

where ^pi is a continuous function on [l,oo), such that 

/oo pOO 
V'i(A)dA = l, J X''^i{X)d\ = 0, k>0, |V'i(A)| < Civ(l + A)-^, N>0 

The existence of such a function was proved in [S] where the extension operator was also introduced. 

In [S] it was proven that this operator is continuous on the Sobolev spaces but it was not proven there 
that it is continuous on the Holder spaces so we must prove this. As pointed out above, we only need 
to prove that it is Holder continuous with respect to the radial variable. 
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First we note that if / G C*' then the extension is in C^. In fact 



(17.16) 



di8xtrU){r) 



fU)(^r - 2A(r - 1))(1 - 2Xy V'i(A) dX, r > 1 



Prom the continuity of d^J and (17.14)-(17.15) it follows that liiRr^i d^. Sxti{f){r) = 5^/(1), that 
£xti{f) is in , and that for A; integer 



(17.17) 



sup|a,^fxii(/)(r)| < Cfcsup|/«(r)| 



Suppose now that A; < a < A; + 1 where k is an integer. We will prove that 

(17.18) sup , c„ sup l-^'"'^' -/'I'MI 

If r < 1 and p <l there is nothing to prove. Also ifr<l<porp<l<r, then Ir — p| > |1 — p| and 
\r — p\>\^ — r\ so in this case, we can reduce it to two estimates with either r = 1 or p = 1. Also it is 
symmetric in r and p so it only remains to prove the assertion when r > p>l. Then we have 

/oo 
[f^^^ (r - 2A(r - 1)) - f^^^ {p - 2\{p - I))) {1 - 2\f ^P^^X) d\ 



< sup ■ 



I rj^f pf I Qi k 



\r- p\ 



a — k 



|(l-2A)'^Vi(A)|dA 



and using the last estimate in (17.15), (17.18) follows. □ 



18. The Nash Moser Iteration. 



At this point, given the results stated in sections 11-14, the problem is now reduced to a com- 
pletely standard application of the Nash-Moscr technique. One can just follow the steps of the proof of 
[AG,H1,H2,K1] replacing their norms with our norms. The main difference is that we have a boundary, 
but we have constructed smoothing operators that satisfy the required properties for the case with 
a boundary. Furthermore, we avoid doing smoothing in the time direction, a similar approach was 
followed in [K2]. Alternatively, one could follow the approach of [Ha], where it is proven that C°° of a 
compact manifold with a boundary is also a tame space, just one small detail is missing which is that 
the the set [0,T] x is not smooth at {0} x and again we get back to the situation were it is 
preferable just to do smoothing in the space directions only. 

We will follow the formulation from [AG] which however is similar to [H1,H2]. The theorem there 
is stated in terms of Holder norms, with a slightly different definition of the Holder norms for integer 
values. However, the only properties that are used of the norms are the smoothing properties, (17.5)- 
(17.8) and the interpolation property (17.3) which we proved with the usual definition, i.e. the one used 
in [HI]. 

Let us also change notation and call $(ii) in (2.13) $(«). Let 

(18.1) |||'"|||a,fe = sup \\u{t,-)\\a,oo + ■■■ + \\DtU{t,-)\\a,oo, || I""!! |a = || k|| |a,0, 

G<t<T 

where ||n(t, ■)\\a are the Holder norms, see (17.1). Proposition 15.1 and Proposition 16.1 now says that 
the conditions and ("^2) below hold: 
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(7^1 ): is twice difFerentiable and satisfies 

(18.2) |p"(ti)(t;i,?;2)|||a < C'„(^|||?;i|||a+^,2|||^^2|IU,2 + |||^^l|IU,2|||^'2|||a+M,2) 

+ Ca\\\ \a+iJ,,2 II |/i,2 |||'J^2 II \n,2, 

where /x = 5, for u,vuV2 G C°°{[0,T],C"°{n)), if 

(18.3) ||k|IU,2<l, Ai = 5 

(W2): Hue C°°{[0,T],C'^m) satisfies (18.3) tlien tiiere is a linear map ij{u) from 
C°°{[0,T],C°°{n)) to C°°([0,r],C°°(n)) such that $'(«)V'(ti) = Id and 

(18.4) IIIV'(^i)5ll|a,2 < Ca{\\\g\\\a+X + MWx ||kl||a+d,2), 

where A = [n/2] + 3 and d = [n/2] + 7. 

Proposition 18.1. Suppose that $ satisfies (Hi), (7i2) and <I>(0) = 0. Let a > iJ,, a > d, a > X + 2//, 
a ^ N. Then 

i) There is neighborhood Wg = {f e C°°([0, T], C°°(n)); |||/|||a+A < S^}, S > 0, such that, for 
f G Ws, the equation 

(18.5) $(n) = / 
has a solution u = u{f) G C^([0, T], C°°(ri)). Furthermore, 

(18.2) ||K/)||U,2<C|||/||U+A, a<a 

In the proof, we construct a sequence Uj G C°°{[0,T],C°°{U)) converging to u, that satisfy 
111^111^,2 < 1 and II IS'iUi 111^^2 < 1; for all j, where Si is the smoothing operator in (18.7). The es- 
timates (18.2) and (18.4) will only be used for convex combinations of these and hence within the 
domain (18.3) for which these estimates hold. 

Following [H1,H2,AG,K1,K2] we set 

(18.7) Ui+i =Ui + dui, 6ui = ip{SiUi)gi, uq = 0, Si = Sq., 0i = Qq1\ 6*0 > 1 
and gi are to be defined so that Ui formally converges to a solution. We have 

(18.8) $(ui+i) - ^{ui) = <^'{ui){ui+i - Ui) + e^' = $'(n,)^(5iU,)5, + e'^ 

= {^'{Ui) - ^'{SiUi))i^{SiUi)gi +gi + e'l = + e'l + gi 

where 



(18.9) e'i = {^'{ui) - ^'{S,u,))5ui 

(18.10) e'l = $(tx,+i) - ^{Ui) - ^'{Ui)5ui 

(18.11) ei = e'i + e!l 
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Therefore 

(18.12) Hui+i) - ^Ui) =ei + gi 

and adding, we get 



i-l 



(18.13) ^Ui+i) = Y,9o+ SiEi + ei + {I- Si)Ei, Ei = J2 

j=0 j=0 

To ensure that $(«,) ^ / we must have 

i 

(18.14) J2gj+SiEi = Sif 

Thus 

(18.15) 9o = Sof, gi = {Si - Si-i){f - Ei_i) - SiCi-i 
and 

(18.16) ^Ui) = Sif + e^ + {I-S,)E, 

Given uq,ui, ...,Ui these determine 5uq, Sui, Sui which by (18.9)-(18.10) determine ei, ej_i, which 
by (18.15) determine gi. The new term Ui+i is the determined by (18.7). 

Lemma 18.2. Assume that |||'Uj|||;:i,2 < 1; || l^ti+i || U,2 < 1 and |||>S'j'Uj|||^^2 < 1- Then 

(18.17) \\\e'i\\\a < Ca(^\\\{I - Si)Ui\\\a+^,,2\\\SUi\^^^^^ 

iUi\\\a+ij.,2\\\{I — Si)Ui\\\^^2\\\SUi\\\^^2 

and 

(18.18) |||e^'|||a<a(|||<5 
Proof. The proof of (18.17) makes use of 

(18.19) {<^'iui) - <P'iSiUi))6ui = [ ^"{SiUi + s{I - Si)ui)iui - Sim, Sui) ds 

Jo 

together with (18.2). Note that from the third term in (18.2) we get a term that is not present in (18.17) 
since it can be bounded by the others using the assumptions. In fact, since |||«i|||;i,2 + |||'S'i?ii|||/ii,2 < 2, 
111(7- S'i)'Ui|||a+^,2|||(i"- 'S'i)ni|||;^,2|||(5'Ui|||;u,2 < 2|||(7- S'i)ni|||„+^,2|||'5ni|||^,2- (18.18) makes use of 

(18.20) Hui+i) - $(ui) - <P'{ui)dui = / (1 - s)<^"{ui + sSui){6ui,Sui) ds 

Jo 

together with (18.2). Here we used that |||'5ui|||a+^,2|||'5^ii|||^,2 < 2|||(5ni|||a+;^x,2|||'^^i^||U,2 □ 

Let a> a and a — fi> 2{a — jj). Throughout the proof Ca will stand for constants that depend on 
a but is independent of n in (18.21). 

Our inductive assumption (i7n) is, 

(18.21) \\\5ui\\\a,2 < 56"^-'' , 0<a<a, i<n 

If n = then if a < a, we have |||5tio|||a,2 < C'alll/IIU+A < C'a^^, so it follows that (18.21) hold for 
n = if we choose 5 so small that CgJ < 9q~". We must now prove that (77„) implies (i7„+i) if 
< 1, where is some constant that only depends on a but is independent of n. 
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Lemma 18.3. // (18.21) hold then fori<n 

i 

(18.22) ^|||5«j|||a,2 <C„(5(min(z,l/|a-a|) + l)(er" + l)> 0<a<a 

3=0 

Proof. Using (18.21) we get E5=o IH'^^.jllU-a < C^a^Ej-o 2^^"""^ ^nd noting that Ej-o^"*^' < 
C(min (1 + l/s,z) + 1), if s > 0, (18.22) follows. □ 

Lemma 18.4. // (i?„), i.e. (18.21), hold and q.> a, then for i <n + l we have 

(18.23) |||ni||U,2 < C„(5(min(z,l/|a-a|) + l)(0r" + l)' 0<a<a 

(18.24) \\\SiUi\\\a,2<CJ{mhl{^^/\a-a\) + lW^-'' + 1), a>0 

(18.25) |||(/- .Si)^illla,2 < C„,5^r"' 0<a<a 

Proof. The proof of (18.23) is just summing up the series Ui+i = Ej=o '^'"i' using Lemma 18.3. (18.24) 
follows from (18.22) using (17.5) for a < a and (17.6) with 6 = d for a > a. (18.25) follows from (17.7) 
with b = a and (18.23) with a = a. □ 

Since we have assumed that a > /U, we note that in particular, it follows that 

(18.26) |||^ii||U,2 < 1 and |||S'itii|||^,2 < 1, for i<n+l if C^5<1. 
As a consequence of Lemma 18.4 and Lemma 18.2 we get 

Lemma 18.5. // (i?n) is satisfied and a > /i, then for i < n, 

(18.27) \\H\\\a < CJ^Or^^"~''\ 0<a<a-n 

(18.28) |||e-'||U < C„(5X"^^""''\ < a < a - /x 

As a consequence of Lemma 18.5 and (17.8) we get 
Lemma 18.6. // (Hn) is satisfied, then for i < n + 1, 

(18.29) III Siei-ilWa < CjHr''^''-^\ a > 

(18.30) \\\{Si - -Si_i)/||U < Caer^\\\f\\\p, a > 

(18.31) 111(7 - S,)/|||„ < CaOr^WlfWy, 0<a<(3 
Furthermore, if a — fx > 2{a — fx): 

(18.32) 111(5, - 5,_i)i^,_i||U < Ca5^er''^"-''\ a > 

(18.33) \\\{I-Si)Ei\\\a<Ca6^e.-^^''-^\ 0<a<a-fi 

Proof. (18.29) follows from (18.27); For a < d — ^ we use (17.5) with b = a and for a > d — /x, wc use 
(17.6) with b = a- n. (18.30) follows from (17.8) and (18.31) follows from (17.7). Now, Ei = EJIo ej 

so by Lemma 18.5 |||£'i|||5,_^ < Ca6'^ E;=o < C^5^6'f~^~^''"~'^'', since we assumed that the 

exponent is positive. (18.32) follows from this and (17.8) with b = a — jj, and similarly (18.33) follows 
from (17.7) with b = a — ji. □ 

It follows that: 
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Lemma 18.7. // (-ffn) is satisfied, a — fi > 2{a — fi), and a> ji then for i <n + l, 

(18.34) < CaSHr^^''-^^ + Caer^wuwy, a > 0. 

Using this lemma and (18.4) we get 
Lemma 18.8. // (i?n) holds, a — ^Ji> 2{a — fi), a > fi, a > d then, for i <n + 1, we have 

(18.35) |||^^i|||a,2 < C„<52e^^-'("-'^^ a > 0. 

Proof Using (18.7), (18.4), (18.34) and (18.24) we get 

(18.36) \\\6u,\\\^,2 < C„(<520f+^-^(«-^) + lll/III^C^"-^ 

+ Caid'et'^'^-^^ + \\\f\\\pe^-^)6{min{z, l/\a-a-d\) + l)(0f+''-" + 1) 

The lemma follows from using that min (i,l/\a — a — d\) + l< C0«/(^,"+'^-" + 1), where C is a constant 
depending on a — d > but independent of i. □ 

If, we now pick /? = a + A, and use the assumptions that A + a < 2[a — /i), and |||/|||a+A < S"^, we 
get that for z < n + 1, 

(18.37) \\\Sui\\\a,2<Ca5^e^-^, a>0. 
If we pick (5 > so small that 

(18.38) CaS < 1, 
the assumption (Hn+i) is proven. 

The convergence of the is an immediate consequence of Lemma 18.2: 

oo 

(18.39) ^ \\\ui+i - Ui\\\a,2 < CaS, a<a 

i=0 

It follows from Lemma 18.6 that 

(18.40) \m'^i)-f\\\a<cj^0r°"^ 

which tends to 0, as z — oo, if a < a + A. 

It remains to prove u G C^([0, T], (il)). Note that in Lemma 18.8 we proved a better estimate 
than {Hn). In fact if we let 7 = 2{a — //) — (a + A) > and a' = a + 7, then || |/|| |a'+A < C implies that 

(18.41) \\\Sui\\\a,2 < CaO^-''' , a>0 

Using this new estimate, in place of {Hn), we can go back to Lemma 18.3-Lemma 18.8 and replace a 
by a' and 5 by 1. Then it follows from Lemma 18.8 that 

(18.42) \\\Sv.\\\a,2 < + aC^^-^lll/lll^ 

and if we now pick 7' = 2{a' — fi) — {X — a') = 27 and a" = a' = a + 27, and use that |||/|||Q,'+-y' < C 
we see that 

(18.43) \\\Sui\\\a,2 < Cjr"\ «>0 

Since the gain 7 > is constant, repeating this process yields that (18.41) holds for any a' and 
hence that (18.39)-(18.40) hold for any a > 0, (with S replaced by 1). It follows that uj is a Cauchy 
sequence in C2([0,T], C'=(0)), for any k, and hence that Uj ^ u e C^{[0,T],C°°{Tl)) and <^{uj) 
f G C{[0,T],C°°m). (18.6) follows from (18.37) with 6^ = |||/|||a+A. This concludes the proof of 
Proposition 18.1. 
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